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scattering rate as flow parameter, we derive formally exact renormalization group flow equations for 
the non-equilibrium self-energies in the Keldysh basis. A simple perturbative truncation of these 
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not suffer from secular terms and gives accurate results even for long times. We demonstrate this 
explicitly within a simple exactly solvable toy model describing a quartic oscillator with off-diagonal 
pairing terms. 
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I. INTRODUCTION 

Quantum mechanical many-body systems out of equi- 
librium pose extraordinary challenges to theory. Al- 
though powerful field-theoretical methods to formulate 
non-equilibrium problems in terms of Green functions 
and Feynman diagrams are available,—^— in practice new 
concepts and approximation strategies are needed. Be- 
cause systems under the influence of external time- 
dependent fields are not time-translationally invariant, 
one has to formulate theoretical descriptions in the time 
domain. Moreover, even if after a sufficiently long time 
the system has approached a stationary non-equilibrium 
state, such a state can exhibit properties which are rather 
different from a thermal equilibrium state. For exam- 
ple, the Fourier transform rik of the distribution function 
at such a non-thermal fixed point can exhibit a scaling 
behavior as function of momentum k which is charac- 
terized by a completely different exponent than under 
equilibrium conditionsj^ii^ In this case a simple pertur- 
bative approach based on the quantum Boltzmann equa- 
tion with collision integrals calculated in lowest order 
Born approximation is not sufficient. The scaling behav- 
ior close to non-thermal fixed points in simple models 
has been studied within a next-to- leading order 
approximation.— This method has also been used to 
study the real-time dynamics of quantum many-body 
systems far from equilibrium^^ 

Another useful method to investigate strongly corre- 
lated many-body systems is the renormalization group 
(RG), which has been extensively used to study the scal- 
ing properties of systems in the vicinity of continuous 
phase transitions^^ Although there are many successful 
applications of RG methods to systems in thermal equi- 
librium, there are only a few examples where RG meth- 
ods have been used to study quantum mechanical many- 
body systems out of equilibrium. Some author o^°'^^ have 
focused on stationary non-equilibrium states, where the 
system is time-translationally invariant so that the quan- 



tum dynamic equations can be formulated in the fre- 
quency domain. On the other hand, the more diffi- 
cult problem of obtaining the time-evolution of quantum 
many-body systems out of equilibrium has been studied 
using various implementations of the RG idea, such as the 
numerical renormalization group approach^ the density 
matrix renormalization groupfi^ real-time RG formula- 
tions in Liouville space?^ and a flow equation approach 
employing continuous unitary transformations J^^i^ In re- 
cent years a number of authors have also applied func- 
tional renormalization group (FRG) methods to study 
many-body systems out of cquilibrium i^^i^^" — While the 
non-equilibrium FRG approach to quantum dots^SriS^ 
has so far only been applied to study stationary non- 
equilibrium states, Gasenzer and Pawlowski22, have writ- 
ten down a formally exact hierarchy of FRG flow equa- 
tions describing the time evolution of the one-particle 
irreducible vertices of interacting bosons out of equilib- 
rium. Using a sharp time cutoff as RG flow parameter, 
they showed that a simple truncation of the FRG vertex 
expansion at the level of the four-point vertex reproduces 
the results of the next-to-leading order l/7V-expansioni^ 

The speciflc choice of the cutoff procedure is a cru- 
cial point when performing RG calculations and differ- 
ent schemes had been proposed in earlier attempts to 
describe non-equilibrium problemsJ^^i^ In this work 
we shall propose an alternative version of the non- 
equilibrium FRG which uses an external out-scattering 
rate as flow parameter. Such a cutoff scheme is closely 
related to the "hybridization cutoff scheme" recently pro- 
posed by Jakobs, Pletyukhov and Schoeller<22i2^i2i Tech- 
nically, we implement this cutoff scheme by replacing the 
inflnitesimal -izirj defining the boundary condition of the 
retarded and advanced Green functions by finite quanti- 
ties ±iA. Given this cutoff procedure, a simple substi- 
tution in the usual hierarchy of FRG fiow equations for 
the irreducible vertices^^ yields the FRG flow equations 
describing the evolution of the irreducible vertices as the 
flow parameter A is reduced. An important property 
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of our cutoff scheme is that it preserves the triangular 
structure in the Keldysh basis and hence does not vio- 
late causality. 

To be specific, we shall develop our formalism for the 
following time-dependent interacting boson Hamiltonian, 



n{t) 
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- y 

2V ^ 



Ik 



5ki+k2+k3+ki,0 

kik2k3k4 

X U{ki, k2; ks, k4)a''_^^^a'_,^^ak^ak^, 



(1.1) 



where creates a boson with (crystal) momentum k 
and energy dispersion e^, the volume of the system is 
denoted by V, and C/(fci, ^2; ^3, ^4) is some momentum- 
dependent interaction function; the minus signs in front 
of the momentum labels of the creation operators in the 
last line of Eq. (|1.1[) are introduced for later convenience. 
As explained in Refs. [2^,|3^ the Hamiltonian (|l.ip de- 
scribes the non-equilibrium dynamics of magnons in or- 
dered dipolar ferromagnets such as yttrium-iron garneli^ 
subject to an external harmonically oscillating microwave 
field with frequency luq. The energy scale 7^ is then pro- 
portional to the amplitude of the microwave field. The 
non-equilibrium dynamics generated by the Hamiltonian 
(jl.ip is very rich and exhibits the phenomenon of para- 
metric resonance for sufficiently strong pumping i^^i^" 

In practice it is often useful to remove the explicit 
time dependence from the Hamiltonian H{t) in Eq. (|l.ip 
by going to the rotating reference frame. The effective 
time-independent Hamiltonian H in the rotating refer- 
ence frame (denoted by a tilde) is obtained as follows: 
Given the unitary time evolution operator U (t) defined 

by 



idtUit) = H{t)U{t), 
we make the factorization ansatz 



with 



(1.2) 



(1.3) 



(1.4) 

where (pk is the phase of 7^ — |7fc|e*'^'=. The time evo- 
lution operator U (t) in the rotating reference frame then 
satisfies 



idtU{t)='HU{t), 



(1.5) 



where the transformed Hamiltonian % does not explicitly 
depend on time, 



H =ul{t)[n{t) ~ idt]UQ{t) 

k ^ 



t \lk\ / t t 

eka'^ak + [^kO-Lk + a,-kak 



5ki+k2+k3+kifi 

kik2k3k4 

xU{ki,k2; k-s, ki)a^_^^a}_^f^^ak3aki, (1.6) 



with the shifted energy 
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(1.7) 



The solution of Eq. (fL5|) is simply U{t) = e"*'^*, so that 
the total time evolution operator of our system can be 
explicitly written down, 



U{t) 



(1.8) 



Throughout this article we shall work in the rotating ref- 
erence frame were the effective Hamiltonian (|1.6p is time 
independent. To simplify the notation we shall rename 
Efc — > Cfc and give all Green functions and distribution 
functions in the rotating reference frame. Explicit pre- 
scriptions to relate these functions in the original and 
the rotating reference frame are given in appendix A. We 
emphasize that the general ERG formalism developed in 
this work remains also valid if the Hamiltonian depends 
explicitly on time. 

The rest of this paper is organized as follows: In Sec. HI] 
we shall define various types of non-equilibrium Green 
functions and represent them in terms of functional inte- 
grals involving a properly symmetrized Keldysh action. 
Due to the off-diagonal terms in our Hamiltonian (|l.ip 
the quantum dynamics is also characterized by anoma- 
lous Green functions involving the simultaneous creation 
and annihilation of two bosons. To keep track of these 
off-diagonal correlations together with the usual single- 
particle correlations, we introduce in Sec. |ll] a compact 
matrix notation. In Sec. lIIII we then derive several equiv- 
alent quantum kinetic equations for the diagonal and off- 
diagonal distribution functions. In Sec. lIVI we write down 
formally exact ERG flow equations for the self-energies 
which appear in the collision integrals of the quantum 
kinetic equations discussed in Sec. IIIII We also discuss 
several cutoff schemes. In Sec. |V] wc then use our non- 
equilibrium ERG flow equations to study the time evo- 
lution of a simple exactly solvable toy model which is 
obtained by retaining in our Hamiltonian (jl.6p only a 
single momentum mode. We show that a rather simple 
truncation of the ERG flow equations yields already quite 
accurate results for the time evolution. Einally, in Sec. lVIl 
we summarize our results and discuss some open prob- 
lems. There are four appendices with additional technical 
details. 



II. NON-EQUILIBRIUM GREEN FUNCTIONS 

Our goal is to develop methods to calculate the time 
evolution of the diagonal and off-diagonal distribution 
functions 



nk{t) = {al{t)akit)) , 
Pk{t) = (a_fc(i)afc(<)). 



(2.1a) 
(2.1b) 



where all operators arc in the Hcisenberg picture and 
the expectation values are with respect to some initial 
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density matrix /5(to) specified at some time t^, 

(...) =Tr[/}(to)...]. (2.2) 

Note that in our Hamiltonian (jl.6p the combinations 
^fc^-fc and a-fcflfc do not conserve particle number, so 
that we should also consider the anomalous distribution 
function (|2.1b[) and its complex conjugate. Our final goal 
is to derive renormalization group equations for the self- 
energy appearing in the collision integrals of the quantum 
kinetic equations for these distribution functions. In or- 
der to do this, it is useful to collect the various types of 
non-equilibrium Green functions into a symmetric ma- 
trix, as described in the following section. 



A. Keldysh (RAK) basis 

In the Keldysh technique^ one doubles the degrees of 
freedom to distinguish between forward and backward 
propagation in time. As a consequence, all quantities 
carry extra indices p £ {-|-, — } which label the branches 
of the Keldysh contour associated with the forward and 
backward propagation. The single-particle Green func- 
tion is then a 2 x 2 matrix in Keldysh space. Alternative 
formulations of the Keldysh technique in combination 
with diagonal and off-diagonal terms are known from the 
theory of non-equilibrium superconductivity, see e.g. js^ . 
In order to formulate the Keldysh technique in terms of 
functional integrals^ it is convenient to work in a basis 
where causality is manifestly implemented via the vanish- 
ing of the lower diagonal element of the Green function 
matrices in Keldysh space. The other matrix elements 
can then be identified with the usual retarded (i?), ad- 
vanced (A) and Keldysh {K) Green functions. Keeping 
in mind that our model has also anomalous correlations, 
we define the following non-equilibrium Green functions, 

9§it,t') = -ie{t-t'){[akit),al{t')]), (2.3a) 
3^(t,t') = ^e(t'-t)(K(t),at(t')]), (2.3b) 
g^{t,t')^-t{{ak{t),al{t')}), (2.3c) 

where [ , ] denotes the commutator and { , } is the 
anticommutator. The corresponding off-diagonal Green 
functions are 

p^it,t') - -te(t-t'){[akit),a^kit')]), (2.4a) 
pt{t,t') = iQ{t' - i)(K(i), a^k{t')]), (2.4b) 
p^(i,t')--«(KW,«-fc(0})- (2.4c) 

In Sec. Ill CI we represent these Green functions as func- 
tional integrals involving a suitably defined Keldysh ac- 
tion. To write the Gaussian part of this action in a com- 
pact form, it is convenient to introduce infinite matrices 
and p^ in the momentum and time labels (where 
X = R,A, K labels the three types of Green functions) 
whose matrix elements are related to the Green functions 



(|2:3aH2:4c)) as follows^ 

[g^]kt,k't' =5k,-k'gkit,t'), (2-5a) 
[p^lfct.fc'f =(5fe,-fc'Pfc (i,i')- (2.5b) 

We have assumed spatial homogeneity so that the Green 
functions are diagonal matrices in momentum space. 
From the above definitions it is easy to show that the 
normal blocks satisfy the usual relations^ 

{g^'V^g^ , (.9^)t = -ff^, (2.6) 

while the pairing blocks have the symmetries 

{p^f=p^ , {p^f=p^. (2.7) 

For each type of Green function, we collect the normal 
and anomalous components into larger matrices, 

A/? „ ( Gaa G^a \ ^ ( £1^ \ fr, q„\ 

'A f G^- \ ( p^ g^ \ 

^[gIgI)-[ {g^r iP^r ) ' 

riK _ f G^a G^a \ ^ ( P^ 9^ \ (0 Rn\ 

^ ^[gl Gf. ) " I -m* -iP'^r ) ' ^'-''^ 

whose blocks G^L/ arc infinite matrices in the momentum 
and time labels. The subscripts cr, cr' e {a, a} indicate 
whether the associated Green functions involve annihila- 
tion operators a or creation operators . We shall refer 
to these subscripts as flavor labels. The symmetries (|2.6p 
and (P?7)) imply 

(G'«)^ = G^ , (G^)^ = G^. (2.9) 

Finally, we collect the blocks (|2.8aff2?8c| into an even 
larger matrix Green function, 

^-[[Gf^ )-[g^ )' ^^-^"^ 

where the superscripts C and Q anticipate that in the 
functional integral approach we shall identify the corre- 
sponding block with correlation functions involving the 
classical (C) and quantum (Q) component of the field, 
see Eqs. (|2.35al2.35dp below. The symmetries im- 
ply that the infinite matrix G is symmetric, 

G = G^. (2.11) 

As emphasized by Vasilie\«^ (see also Refs. [H,!!!!) the 
symmetrization of the Green function and all vertices 
greatly facilitates the derivation of the proper combina- 
torial factors in perturbation theory and in the functional 
renormalization group equations. The definitions (|2.5ai 
I2.5bl I2.8c[) imply that at equal times and vanishing total 
momentum the matrix elements of the Keldysh block G^^ 
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contain the diagonal and off-diagonal distribution func- 
tions defined in Eqs. (|2.1al2.1bp . 



Pkit,t) g^it,t) 

Pk(t) nk(t) + \ 



-2i 



(2.12) 



For later reference we note that the inverse of the matrix 
G in Eq. (|2.10p has the block structure 



G 



(2.13) 



where the products in the lower diagonal block denote 
the usual matrix multiplication, i.e., 

o"l fci 

(2.14) 

The symmetry relations p.9p guarantee that the lower 
diagonal block in Eq. (|2.13p is again symmetric. To in- 
troduce a matrix F in flavor space which contains both 
the normal and the anomalous components of the distri- 
bution function, we parametrize the Keldysh block in the 
form 

qK ^QRp^ _ ^pTQA^ (2.15) 
where the antisymmetric matrix Z = —Z^ is defined by 



Z = Z®1 



1 

-i 



(2.16) 



Here Z is the following antisymmetric 2 x 2-matrix in 
flavor space. 



Z = i(72 



1 
-10/' 



(2.17) 



and 1 is the unit matrix in the momentum- and time- 
labels, i.e. [l]fct.fc't' = Sk,k'S{t — t'). Substituting 
Eq. (|2.15p into (|2.13p . the lower diagonal block of the 
matrix (|2.13p can be written as 

[G-i]W = -(G«)-iG^(G^)-i 

= {G^y^ZF^ - FZ{G^)-\ (2.18) 

so that Eq. (|2.13p takes the form 

_ f (G^^-i 



G 



V (G^)-i (G«)-iZF^ - FZ{G^y 



(2.19) 

In Sec. IIIDI we shall explicitly calculate the matrix ele- 
ments of F in the non-interacting limit, see Eqs. (|2.68[ 
[270)1 below. 

Introducing the self-energy matrix S in all indices via 
the matrix Dyson equation 



G-^ = G„ 1 



(2.20) 



the non-equilibrium self-energy in the Keldysh basis ac- 
quires the form 



S = 





[SI fSl 





( ' 











(2.21) 



The sub-blocks contain the normal and anomalous self- 
energy matrices 



yR 
■^aa aa 
^R vfl 



yH yR_ 
aa aa 



'-^A yA_ 

-'aa aa 
'-^A s^A 



yA y/l_ 
aa aa 



yK 
-'aa aa 



yK y« 

aa aa 



{a^Y (TT^) 
-{a'<Y H^") 



(2.22a) 
(2.22b) 
(2.22c) 



The Dyson equation ()2.20p and the symmetries 
(|2.6l2.7l2.9p imply that the sub-blocks satisfy the sym- 
metries 



(a«)t=<T- , (a-) 



and 



R\T 



(tt-) 



K\T 



(^-) 



(2.23) 
(2.24) 



Since the self-energy blocks satisfy the same symmetry 
relations as the Green functions p.6l2.7p , the symmetries 
p.9p hold also for the self-energy blocks in Keldysh space. 



{t"f = t^ , (E^)^ = E^. (2.25) 

The full self-energy matrix is therefore again symmetric, 

S = YF. (2.26) 

In the presence of interactions the lower-diagonal block 
of the inverse propagator is given by the negative of the 
Keldysh component of the self-energy. 



(2.27) 



B. Contour basis 



In the Keldysh technique all operators are considered 
as functions of the time-argument on the Keldysh con- 
tour. The time contour runs in real-time direction from 
some initial time to some upper limit t>, which is 
larger than all other times of interest and which is slightly 
shifted in the upper positive imaginary branch of the con- 
tour, and then back to to in the lower, negative imag- 
inary branch. Alternatively, all time integrations can 
be restricted to the interval [to)i>] and one can keep 
track of the two branches of the Keldysh contour us- 
ing an extra label p = ±, where p = -\~ corresponds to 
the forward part of the contour and p = — denotes its 
backward part. In the functional integral formulation 
of the Keldysh technique^ the bosonic annihilation and 
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creation operators are then represented by pairs of com- 
plex fields ak p{t) and ak,p{t) carrying the contour label 
p. The contour ordered expectation values of these fields 
define four different propagators G^p , which are related 
to the usual time-ordered (G^), anti-time-ordered (G^), 
lesser (G^) and greater (G^) Green functions and their 
RAK-countcrparts as foUows^i^ 



Rig. = m. 



G< \ _ f G++ G+- 
G> G'^ ) ^\ G-+ G— , 

(2.28) 

where the transformation matrix R has the block struc- 
ture 



V2 U -I 



(2.29) 



Here / is the unit matrix in the flavor-, momentum-, and 
time labels, i. e., [I]a kt,a'k't' = Sa,cr'Sk,k'S{t-t'). The ma- 
trix equation p.28p implies for the blocks in the contour 
basis 



GPP 



p'G^+pG' 



G 



K 



(2.30) 



where p,p' e {+,—}. From Eq. (|2.30p one easily verifies 
the inverse relations. 



pp' 
pp' 

^Gfc^^Y^OPP'^ 



= J2pp'GPp'. 

pp' 



(2.3fa) 
(2.3fb) 
(2.3fc) 
(2.3fd) 



The corresponding relations for the self-energy are 



which gives 
Y,pp' 



pt^ +p'T.'^ +pp'Y. 



R 







and the inverse relations 



pp' 

pp' 



pp' 



= ^EPP'. 
pp' 



R, 
(2.32) 

(2.33) 

(2.34a) 

(2.34b) 

(2.34c) 
(2.34d) 



Functional integral representation of Green 
functions in the RAK-basis 



To define the proper boundary conditions in the func- 
tional integral formulation of the Keldysh technique,— it 
is convenient to work in the RAK basis. The transition 
from the contour basis to the RAK basis is achieved by 
introducing the classical (G) and quantum (Q) compo- 
nents of the fields, which are related to the corresponding 
fields flfc -t and afc.± in the contour basis via 





f 


[dk,- 


f(0 


+ flfc - 




(2.35a) 


-4{t) = 


1 


[ak,- 




+ flfc,- 




(2.35b) 


= 


f 

71 


[dk,- 


f(0 


- flfc,- 




(2.35c) 


-4(t) = 


1 

71 


[flfc,- 




- flfc - 




(2.35d) 



Introducing a four-component "super-field" 



<i>Q{k,t) 



(2.36) 



the matrix elements of the symmetrized matrix Green 
function G defined in Sec. Ill Al can be represented as the 
following functional average. 



z[G]^^;,,,;,,,, ^iG^^:(fci,fcY) = ($^(fe<)ci>^;(fc't')) 

2?[$]e''5[*l$^(M)$^;(fcY). (2.37) 



For simplicity, we shall assume throughout this work that 
the initial state at t = to is not correlated. In this case the 
corresponding Gaussian initial density matrix docs not 
explicitly appear in the above Keldysh action S\^\^ but 
enters the dynamics via the initial conditions for the in- 
dependent one- and two-point functions.— i^^— In prin- 
ciple, non-Gaussian initial correlations can be considered 
either explicitly by assuming a non-Gaussian initial den- 
sity matrix, or implicitly by introducing multiple sources 
in the generating functional for the n-point functions]^ 
The latter approach will also change the matrix structure 
of the theory; for instance if the system is initially in ther- 
mal equilibrium the Keldysh matrices expand from a 2 x 2 
to a 3 X 3 form. As shown in Ref. [s^, initial correlations 
can lead to additional damping effects which modify the 
amplitude and phase of the oscillatory evolution at inter- 
mediate times. In the absence of initial correlations the 
form of the Keldysh action in Eq. (|2.37|) can be di- 
rectly obtained from the corresponding Hamiltonian, so 
that it has has the following two contributions. 



^[$] = 5o[*] + ^i[$] 



(2.38) 
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where after proper symmetrization the Gaussian part 
Sq [$] can be written as 

(T(t' W kk' 

X <I>^(fct)[Go-i]^^;.,,fc,,,$^;(fcY). (2.39) 



Here Gq^ is the non-interacting inverse Green func- 
tion matrix in Keldysh space, which is associated with 
the non-interacting part of the Hamiltonian (|1.6p . The 
matrix G,^^ has the same block structure as G~^ in 
Eq. p.l3p . with retarded and advanced blocks given by 



{G^)-'=D-ii^Z , (G^)-' =D + tr^Z, (2.40) 



A\-l 



where the antisymmetric matrix Z is given in Eq. (|2.16p , 
and the symmetric matrix D is defined by 



kt.k't' 



(2.41) 



Here S'{t) ~ ^'^(^) the derivative of the Dirac 6- 
function and is the following matrix in flavor spaced 



Mfc = 



efc |7fc| 
-|7fe| -ffc 



(2.42) 



Recall that we are working in the rotating reference frame 
where we have redefined Cfc = — ^ — )■ efc,see Eq. (|1.7p . 
Keeping in mind that i5'{t—t') = —i5'{t'—t), it is obvious 

that [{G^rY = (Go )"\ in agreement with Eq. 

Although the Keldysh bl ock [ Gp of the inverse 
Gaussian propagator in Eq. (|2.39p vanishes in continuum 
notation, it is actually finite if the path integral is prop- 
erly discretizedii It is, however, more convenient to stick 
with the continuum notation and take the discretization 
effectively into account by adding an infinitesimal regu- 
larization 77. To derive this regularization, we note that 
the relations Gq^^Go = GqGq"^ = I imply that in the 
non-interacting limit the Keldysh block satisfies 



DG{^ = G^D = 0. 



(2.43) 



Introducing the non-interacting distribution matrix Fq as 
in Eq. ()2.15p this implies for ry 0, 



DG^D = FqZD - DZF^ 



0. 



(2.44) 



Using Eq. p.lSp we thus obtain for the lower diagonal 
block of Gg"^ in the non-interacting limit. 



[Go 



= {G^)-'ZFS-FoZ{G^)-' 
= irj{Fo + P;^) = 2ti^Fo, 



(2.45) 



where we have used the fact that Fq is symmetric, which 
is easily verified by explicit calculation, see Sec. Ill D[ The 
lower diagonal block of Gq ^ is thus a pure regularization, 
which guarantees that in the non-interacting limit the 
functional integral (|2.37p is well defined. In the presence 
of interactions the Keldysh component of the self-energy 
is finite due to Eq. (|2.27p . so that in this case the in- 
finitesimal regularization (|2.45p can be omitted. 

To write down the interaction part of the Keldysh 
action associated with the interaction part of the 
Hamiltonian (jl.ip . one should first symmetrize the 
Hamiltonia n"*^^ -^^ before formally replacing the operators 
by complex fields, since a and are treated symmetri- 
cally. Noting that the symmetrized product of n bosonic 
operators ^1, . . . , An is defined by 



{Ai...An} = -^y2Ap,.. 



(2.46) 



where the sum is over all n\ permutations, we have 

+Sk2MWk,akA + Sk2MWk,<^k-A 



5ki,k3^k2,ki + 5k2M^ki,ki ■ (2.47) 



The quadratic terms on the right-hand side of Eq. (|2.47p 
lead to a time-independent first-order shift in the bare 
energy dispersion. 



efc^efc--^^C/(-fc,-fc';fc,fc') 



{2A 



This shift can be absorbed by re-defining the energy in 
the matrix Mk introduced in Eq. (|2.42p . The first term 
on the right-hand side of Eq. (|2.47p leads to the following 
interaction part of the Keldysh action in the RAK-basis, 



S'i[<I>]=-^ ^ J dt Ski+k2+k3+k4,oU{ki,k2;k3,k4^) 

ki ^2^3/24 

X |$c^(fcii)*?(fc2t) [^^(M)*? (M) + <i>?(fe3t)$?(fe4i) 

+ {k,t)^^ik2t) + $?(fcit)$?(fc2i)l *a (fc3t)*?(M) 



(2.49) 



To eliminate complicated combinatorial factors in the ERG flow equations derived in Sec. IIVI it is convenient to 
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symmetrize the interaction vertices in Eq. p.49p with respect to the interchange of any two label a^^i'^^'^^ and write 

(2.50) 



(7i...<74 A1...A4 ki ...k^ 



where the interaction vertex U^^^^^^^^{ki,k2,k3,k4) is 
symmetric with respect to any pair of indices. Up to 
permutations of the indices, the non-zero vertices are 

C/a?af(fcl,fc2,fc3,fc4) 
= f/????(fcl,fe2,fe3,fc4) 

= U^'ia?iki,k2,k3,ki) 

= UiTa'^{k^,k2, fcg, fc4) = U{k^M;k^, ^4). (2.51) 

D. Non-interacting Green functions 

To conclude this section, let us explicitly construct the 
2x2 matrix Green functions in flavor space in the non- 
interacting limit. In general, we define the Green func- 
tions in flavor space in terms of the matrix elements 

[G^]fc^-fct' -G^(fc,t,t') , X = R,A,K. (2.52) 

In the absence of interactions, we can obtain explicit 
expressions for these Green functions. Then the non- 
interacting part of the Hamiltonian (|1.6p reduces to^° 



^o = E 



\lk[ 
2 



O—kO-k 



(2.53) 



The retarded and advanced matrix Green functions in the 
non-interacting limit are now easily obtained. Consider 
first the retarded Green function, which satisfies 

idtG^{k,t,t') ^ S{t - t')Z + MkG^{k,t,t'). (2.54) 

The solution of Eq. (|2.54p with proper boundary condi- 
tion is 



G^(fc, t, t') = -iQ{t - <')e"**^'=^*"*'^^- 
The matrix exponential is 



I cos{pkt) — iMk 



sm{pkt) 



fJ-k 



where / is the 2x2 unit matrix, and 



|7fcp if kfcl > l7fe| 



Wllkl'-e't if |7fc| > kfel 



(2.55) 



(2.56) 



(2.57) 



Using the symmetry relation (G^)^ = given in 
Eq. (|2.9p . we obtain for the corresponding advanced 
Green function. 



G^{k,t,t') = G^{~k,t',t)' 



-iQ{t'-t)Z^e-'^'S<^*'-*\ 
(2.58) 



r 



Using the identities 



Z' = 

z-i = 

Ml = 



-I, 

= -Z, 
ZMkZ = 

_^-iMh{t-t') ^ 



Z^MkZ, 



we may also write 

G^(fc, i, t') = te{t' - Oe-*^^'=^*-*')Z. 



(2.59a) 
(2.59b) 
(2.59c) 
(2.59d) 



(2.60) 



Because the retarded and advanced Green functions de- 
pend only on the time difference, it is useful to perform 
a Fourier transformation to frequency space. 



G^{k,u) = I dte"^'G${k,t,0) 



iujt r~iX I 



(2.61) 



Substituting Eqs. (|2.60p and (|2.55p into this expression 
and representing the step-functions as 



m = 

it is easy to show that 

G^ik^Lo)-- 
G^{k,Lu) = 

or explicitly. 



27ri Lj' — ill 



[uj - Mk + ivr^Z, 
[uj-Mk-ivr^Z 

z[oj+A4-tir\ 



(2.62) 

(2.63a) 
(2.63b) 



Gr(fc,a.) = 



{oj ± iriy 



I7fc| 



-|7fc| u}±ir-i + ek 



(2.64) 



Next, consider the Keldysh component Gq {k,t,t') of 
our 2x2 matrix Green function in flavor space. It satisfies 
the matrix equations 

idtG^ik,t,t') = MkG^{k,t,t'), (2.65a) 
idt'G^{k,tX) G^ {k,t,t')Ml , (2.65b) 

and hence 

{idt+idt>)G^[k, t, t') = MkG^{k, t, t')+G^{k, t, t')Ml. 

(2.66) 

These equations are solved by 



(fc, i, t') = e-'^'<^*G^ (fe, 0, 0)e 



-iMit' 



(2.67) 



8 



with an arbitrary initial matrix Gq {k, 0, 0) which defines 
the distribution functions at t = 0. To expUcitly con- 
struct the distribution matrix Fq defined via Eq. ()2.15p . 
we note that in the non-interacting limit the distribution 
matrix is diagonal in time, 

[Fo]kt,k't' ^ Sk.-k'Fo{k,t,t') = Sk.-k'Sit - t')FQ{k,t), 

(2.68) 

so that Eq. p.lSp reduces to the 2x2 matrix relation 

G^{k, t, t') = G^{k, t, t')Fo{k, t')Z - ZFo{k, t)G^{k, t, t') 
= -ie{t - t')e~'^'''^'-''^ZFo{k, t')Z 

^iQ{t' -t)ZFo{k,t)Ze-'^'^'^''-'\ (2.69) 

It follows that in the non-interacting limit the time- 
diagonal element FQ{k, t) of the distribution matrix con- 
tains the normal and anomalous distribution functions 
defined in Eqs. (|2.1a[ I2.1bp in the following way, 

(2.70) 

Combining this relation with Eq. (|2.66p , we see that our 
diagonal distribution function matrix satisfies the kinetic 
equation 

idtFoik,t) = -MlFo{k,t) - Fo{k,t)Mk. (2.71) 

Note that the non-interacting Keldysh Green function 
(|2.69p can also be written as 

G^{k,t,t') = ~i[G^{k,t,t')ZG^{k,t',t') 

-G^ik,t,t)ZG^ik,t,t% (2.72) 

which relates the matrix elements of the Keldysh Green 
function at different times to the corresponding equal- 
time matrix elements. 

III. QUANTUM KINETIC EQUATIONS 

From the Keldysh component of the Dyson equation 
we obtain quantum kinetic equations for the distribution 
function. In this section we shall derive several equivalent 
versions of these equations. Although matrix general- 
izations of quantum kinetic equations are standard^i^i^ 
we present here a special matrix structure of the kinetic 
equations which takes into account off-diagonal bosonic 
correlations. 

A. Non-equilibrium evolution equations for 
two-time Keldysh Green functions 

In order to derive quantum kinetic equations we start 
with the matrix Dyson equation (|2.20p . which can be 
written as 

(Go 1 - S) G = I. (3.1) 



This "left Dyson equation" is equivalent with the follow- 
ing three equations for the sub-blocks, 

[(G^)-i - ±^]G^' = /, (3.2a) 
[(G^)-i - ±^]G^ = /, (3.2b) 

[(G^)-i-S^]G^ = E^G^, (3.2c) 

where / is again the unit matrix in the flavor-, 
momentum- and time labels. Alternatively we can also 
consider the corresponding "right Dyson equation" , 

G (Go ^ - S) = I, (3.3) 

which implies the following relations. 



G^[(G^)-i - E^i 




(3.4a) 


G^[(G^)-i -S-^] 


= 1 


(3.4b) 


G^[(G^)-i -E^] ^ 


G^E^. 


(3.4c) 



In order to solve the coupled set of equations p.2aH3.2c)) 
and p.4al[5Tlc)) , it is sometimes useful to rewrite them as 
integral equations. Therefore we should take into account 
that in the non-interacting limit the Keldysh self-energy 
is actually an infinitesimal regularization, E^ = —2i'qFo, 
see Eq. (|2.45p . In the non-interacting limit the Keldysh 
component therefore satisfies 

{G§)-'Gl^ = ~2zrjFoGt (3.5) 

which is equivalent with 

G^ = -2if]G^FoG^. (3.6) 

Using the integral forms of the Dyson equations 
p.2al3.2bp for the retarded and advanced Green func- 
tions, 

G^ = G^ + G^E^G^, (3.7a) 
G^ = Gi + G^t'^G^, (3.7b) 

the equation (|3.2cp for the Keldysh-block can alterna- 
tively be written in integral form as 

G^ = G^E«G^ + G«S^G^ 

= G^ + G^t^G^ + G^t^G^ + G^t^G^. (3.8) 

Given an approximate expression for the self-energies, 
the integral equations (|3.7l3.8p can be solved by an ap- 
propriate iteration to obtain the time evolution of the 
Keldysh Green function. Here we follow a different strat- 
egy which is similar to the one proposed in Ref. [ll|. 
We represent the inverse propagators (|2.40p as differen- 
tial operators to derive evolution equations in differential 
form. We will see that the resulting initial value prob- 
lem allows for approximate solutions which manifestly 
preserves causality. Using the advanced and retarded 
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components of the Dyson equations (|3.2al3.2bl3.4al3.4b)) 

and keeping in mind that by translational invariance the 
matrix elements in the rotating reference frame are, 



with X ^ R, A, K , wc obtain, 



(3.9a) 
(3.9b) 



pt/t' 

+ / dtiZJ:^/'^{k,t,ti)G^/'^{k,ti,t'), (3.10a) 

J to 

idt'G^^'^ik, t, t') - G^/^(fe, t, t')Ml = -ZS{t - t') 
ft/t' 

- / dtiG^^^{k,t,ti)i:^^'^{k,ti,t')Z. (3.10b) 

Jtn 



In the same way we obtain from (|3.2cl I3.4cp the following 
kinetic equations for the Keldysh component, 

idtG'^ (k,t,t') ^ MkG^ {k,t,t') 

rt' 



f dhZY.^{k,t,ti)G'^{k,ti,t') 

J to 

[ dhZ^^{k,t,ti)G^{k,h,t'), (3.11) 

Jtn 



and 



idt'G^{k,t,t')-G^{k,t,t')Ml 



[ dtiG^{k,t,ti)^^{k,ti,t')Z 

J to 

f dtiG^{k,t,ti)^^{k,ti,t')Z. (3.12) 

Jtn 



I to 

Finally, in order uniquely define the solution of the 
set of coupled first-order partial differential equations, 
proper boundary conditions for the Green functions have 
to be specified. From the definitions (|2.3al2.3bp and 
(|2.4al2.4bp of the advanced and retarded propagators we 
find that for infinitesimal 77, 



G^{k,t,t-7j) = -iZ, 
G'^{k,t,t + 7]) = iZ, 



(3.13a) 
(3.13b) 



and that the Keldysh Green function should reduce to 
the matrix G'^(fc,0,0) at the reference time t = <' = 0. 
Note that we have not made any approximations so far 
and the time evolution is exact provided we insert the 
exact self-energies. The evolution equations are causal 
by construction, since no quantity in the collision inte- 
grals on the right-hand side depends on future states. 
Interpreting the time derivatives as finite-difference ex- 
pressions, the solution can be obtained by stepwise prop- 
agating the equations in t and f direction. Note that in 
the non-interacting limit where all self-energies and col- 
lision integrals vanish, our kinetic equations p.lll3.12p 



correctly reduce to the equation of motion for the non- 
interacting Keldysh Green function given in Eq. (|2.66p . 

For open systems coupled to an external bath, it is 
sometimes convenient to move some of the terms on the 
right-hand side of the kinetic equations (|3.11l3.12p to the 
left-hand side such that the remaining terms on the right- 
hand side correspond to the "in-scattering" and the "out- 
scattering-rate" in the Boltzmann equation. To achieve 
this, we introduce the average (mean) and the imaginary 
part of the retarded and advanced self-energies,— 

l.^p (3.14) 



The inverse relations are 



2- ' - - ■ 2^'- ^'-''^ 
A similar decomposition is also introduced for the re- 
tarded and advanced Green functions, 

G^^ = i[G'« + G^] , G^ = t[G^-G^], (3.16) 

so that 



G^ = G^' 



2^ 



G^ = G 



M 



2^ ■ 



(3.17) 



Subtracting the Keldysh component of the left and right- 
hand sides of the Dyson equations p.2cl3.4cp . we obtain 
the (subtracted) kinetic equation 

with 



D 



(3.19) 



The collision integrals are represented by symmetric ma- 
trices, 



G'" 

^out 



-[zt^^G' + G^±^Z], 
^[Z±^G^ + G^^±^Z], 



(3.20a) 
(3.20b) 



and correspond to the usual "in-scattering" and "out- 
scattering" term in the Boltzmann equatioui^ The kinetic 
equation (|3.18p generalizes the subtracted kinetic equa- 
tion given in Ref. [Bj to matrix form, which includes also 
off-diagonal correlations. In equilibrium both terms on 
the right-hand side of Eq. p.lSp cancel. 



B. Evolution equations for the equal-time Keldysh 
Green functions 



Keeping in mind that in analogy with Eq. (|2.70p the 
diagonal and off-diagonal distribution functions are con- 
tained in the matrix F{k,t), 



F{k,t) = iZG'^{k,t,t)Z 



-2pl{t) 2nkit) + l\ 
(3.21) 
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we only need to calculate the time evolution of the equal- 
time Keldysh Green function in order to obtain the dis- 
tribution function. Adding Eqs. p.ll|) and p.l2|) and 
using 

{dt + dt')G''{k,t,t%,^^ = dtG''{k,t,t), (3.22) 

we arrive at the evolution equation for the equal-time 
Keldysh Green function, 

idtG^ik, t, t) - MkG^{k, t, t) - G^{k, t, t)Ml 
= [ dti[ZY.^ {k,t,ti)G^{k,ti,t) 

J to 

-G"{k,t,ti)T.^ {k,ti,t)Z] 
+ [ dti[Z^^{k,t,ti)G^\k,ti,t) 

J to 

-G^ik,t,ti)^'^{k,ti,t)Z]. (3.23) 

Although the left-hand side of Eq. p.23p involves the 
Keldysh Green function only at equal times, the inte- 
grals on the right-hand side depend also on the Green 
functions at different times (notice the implicit depen- 
dence via the self-energies), so that Eq. (|3.23p is not a 
closed equation for G^ {k,t,t). In principle, one has to 
solve the more general equations (|3.1ip and (|3.12p which 
fully determines G^ {k,t,t') for all time arguments. An 
alternative strategy is to close the kinetic equation p.23p 
for the equal-time Keldysh Green function by approxi- 
mating all Keldysh Green functions at different times on 
the right-hand side in terms of the corresponding equal- 
time Green function. This is achieved by means of the 
so-called generalized Kadanoff-Baym ansata^ (GKBA) 
which is one of the standard approximations to derive 
kinetic equations for the distribution function from the 
Kadanoff-Baym equations of motion for the two-time 
non-equilibrium Green functions.— i"^^'^^ For bosons with 
diagonal and off-diagonal correlation, the GKBA ansatz 
reads 

G^{k,t,t') « -i[G'^{k,t,t')ZG^{k,t',t') 

-G^{k,t,t)ZG'^{k,t,t')], (3.24) 

which assumes that the exact relation Eq. (|2.72p be- 
tween non-interacting Green functions remains approx- 
imately true also in the presence of interactions. Note 
that Eq. (|3.24p is a non-trivial 2x2 matrix relation in 
flavor space. In appendix B we discuss the approxima- 
tions which are necessary to obtain the GKBA from the 
exact equations of motion (|3.11l3.12p for the two-time 
Keldysh Green function. 

IV. NON-EQUILIBRIUM FUNCTIONAL 
RENORMALIZATION GROUP 

The functional renormalization group (ERG) has been 
quite successful to study strongly interacting systems in 
equilibrium, see Refs. psi . l46j for reviews. In contrast to 



conventional renormalization group methods where only 
a finite number of coupling constants is considered, the 
FRG keeps track of the renormalization group flow of en- 
tire correlation functions which depend on momentum, 
frequency, or time. In principle, it should therefore be 
possible to calculate the non-equilibrium time evolution 
of quantum systems using FRG methods. Recently, sev- 
eral authors have generalized the FRG approach to quan- 
tum systems out of equilibrium] In particular, 
Gasenzer and Pawlowski^S have used FRG methods to 
obtain the non-equilibrium time evolution of bosons. 

Given the Keldysh action defined in 
Eqs. (|2.38I2.39I2.50P it is straightforward to write 
down the formally exact hierarchy of FRG flow equa- 
tions for the one-particle irreducible vertices of the 
non-equilibrium theory, which we shall do in the follow- 
ing subsection. The real challenge is to devise sensible 
cutoff schemes and approximation strategies. We shall 
address this problem in Sec. |V] using a simple exactly 
solvable toy model to check the accuracy of various 
approximations. 

A. Exact FRG flow equations 

To begin with, we consider an arbitrary bosonic many- 
body system whose Gaussian action is determined by 
some inverse matrix propagator G,^^, which we modify 
by introducing some cutoff parameter A, 

Go ^ Go,A. (4.1) 

Depending on the problem at hand, different choices of A 
may be appropriate. For systems in thermal equilibrium 
it is usually convenient to choose A such that it removes 
long-wavelength or low-energy fluctuations.™ In order to 
calculate the time evolution of many-body systems, other 
choices of A are more appropriate. For example, Gasen- 
zer and Pawlowski have proposed that A should be iden- 
tified with a time scale r which cuts off the time evolution 
of correlation functions at long timesi^ In Sec. lIVB^ we 
shall propose an alternative cutoff scheme which uses an 
external "out-scattering rate" as RG cutoff. 

Given the cutoff-dependent Gaussian propagator (|4.ip , 
the generating functional of all correlation functions de- 
pends on the cutoff. By taking the derivative of the gen- 
erating functional Fa [$] of the irreducible vertices with 
respect to the cutoff, we obtain a rather compact closed 
functional equation for rA[<I'] which is sometimes called 
the Wetterich equation4i Formally the Wetterich equa- 
tion is valid also for quantum systems out of equilibrium 
provided we use the proper non-equilibrium field theory 
to describe the system<^ By expanding the generating 
functional rAi*!"] in powers of the fields we obtain the 
one-particle irreducible vertices of our non-equilibrium 
theory, 

rA[*] = E^/ ■ I ril...„„ (4.2) 

„=0 "'"I 
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Here the collective labels ai,a2,... stand for all la- 
bels which are necessary to specify the fields 3^1^ For 
the boson model defined in Sec. HI] the collective label 
a = {X,a,k,t) represents the Keldysh label A S {C,Q}, 
the flavor label a € {a, a}, as well as the momentum 
and time labels k and t. The corresponding integration 
symbol is 



EEE 



dt. 



(4.3) 



The exact FRG flow equations for the irreducible vertices 
can be obtained from the general FRG flow equations 
given in Refs. [l^ll^li^ by making the following substi- 
tutions to take into account the different normalization 
of the action in the Keldysh formalism, 

A,ai...Q„ ~^ A,Qi...C(„' 

Ga -«Ga , Ga —iGa, 
where the single-scale propagator is given by 

Ga = -GaI^aGq a)Ga. 



(4.4) 
(4.5) 

(4.6) 



This definition implies that the blocks of the single-scale 
propagator for general cutoff are 



[Ga 


CQ 


= Gf = 




-']Gl 


(4.7a) 


[Ga 


QC 


— Lta — 


-GiidAiGoA)' 


'']Gl 


(4.7b) 


[Ga] 


CC 


AiK 
— Ua — 


-G^[dA{Go,Ar 







G^[5a(Go^)-']G^ 
G^{dA[G-\]QQ)Gi 



(4.7c) 



If the expectation values of the field components <f>Q van- 
ish in the absence of sources, the exact FRG fiow equation 
for the irreducible self-energy (two-point function) is 



^aF^'' = - 



/3l JP: 



[Ga]/3i/3oF 



(4) 

2^ A,/32/3iaia2 ' 



(4.8) 



while the four-point vertex (effective interaction) satisfies 



^ArA, 



(4) 



[G 



A\fSil32^ A,l32l3iaia2a3C 



(6) 



/3i J 132 J [33 -I p. 
-(4) 



l[3i J 132 
[GA]/3i/32[GA]fe/34 



-(4) 



A,/32&a3a4 A,/34/9iQiQ2 



F 



(4) 



-(4) 



A,02/33aiC(2 A,;34/33a3C(4 



-t-(ai O 02) + («! "i-> O-a) 



(4.9) 



If, in the absence of sources, the field has a finite ex- 
pectation value $° ^ 0, it is convenient to re-define the 
vertices in the functional Taylor series (|4.2p by expanding 
in powers of (5$ = $ - 

rA[<f] = f:;Jy / ••■/ ri'l....J*°)^*.4 • • "^'i'.,.. 

„=0 "'"1 

(4.10) 



Then, the odd vertices (n = 1,3,5,...) are in general 
also finite. The requirement that the one-point vertex 
vanishes identically leads to the following flow equation 
for the fleld expectation value^^i^^i^ 



{dA^l)[Gl%o. + ^l[dAG-XW 



[3i J 13, 



[G 



A\PiP2^A',l32l3ia- 



(3) 



.11) 



Moreover, the FRG fiow equation for the two-point ver- 
tex contains additional terms involving the three-point 
vertex. 



* / / [r 1 

r, / / [*^AJ/3i02^ A,/32^lQlQ2 
^ Jl3l J 02 

I (5A$^)ri%„^„^ 
J 



+^1111 \-^^hl0'^[^A]f32f3;,^A!p^02al'^A.'(3',|3' 

hi J 02 J0'i 



.(3) 



(4.12) 



B. Cutoff schemes 



1. General considerations 

The crucial point is now to identify a sensible flow pa- 
rameter A. Since we are interested in calculating the time 
evolution of the distribution function at long times the 
flow parameter should be chosen such that for sufficiently 
large A the long-time asymptotics is simple. This is the 
case if A is identified with a scattering rate which intro- 
duces some kind of damping. This strategy was already 
implemented by Jakobs et al^ in their recent FRG study 
of stationary non-equilibrium states of the Anderson im- 
purity model. Formally, such a cutoff can be introduced 
by replacing the infinitesimal imaginary part rj appear- 
ing in the retarded and advanced blocks of the inverse 
matrix propagator Gq ^ given in Eqs. p.40p by a finite 
quantity A. This amounts to the following replacement of 
the inverse retarded and advanced propagators by cutoff- 
dependent quantities. 



{G^.a)-' = {G'oT'-^^Z, (4.13a) 



(G^)-' ^ (G^a)"' = (G^r' + *AZ. (4.13b) 

Explicitly, the cutoff-dependent retarded and advanced 
Green functions are then 

G^^^{k,t,t') - Go«(fc,t,i')e-^^*-*'\ (4.14a) 
Gi^{k,t,t') = G^{k,t,t')e-^^''-'\ (4.14b) 

As far as the QQ-component of the inverse free prop- 
agator is concerned (which for infinitesimal 77 is a pure 
regularization) , we set 



(4.15) 
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where the distribution matrix -F^.a will be further 
specified below. Defining the cutoff-dependent non- 
interacting distribution matrix Fo,A as in Eq. (|2.15p . we 
have 



(4.16) 



and 



= DZFo.A - Fo^aZD + 2iAFo,A. (4.17) 



Comparing this with Eq. (|4.15p . we see that our cutoff- 
dependent distribution matrix satisfies 

DZFo,A ~ FoaZD + 2iA(Fo.A - F,,a) = 0. (4.18) 

For our cutoff choice given in Eqs. (|4.13al4.13bp we have 

9a(G?a)^' = -^Z, (4.19a) 



9a(G^a) 



*,Aj 



(4.19b) 
(4.19c) 



so that in this scheme the blocks of the single-scale prop- 
agator (|4.7p are 



— T-'-'A^'-'A^ 
(jA — —^'-^A^L'A: 



G 



K 



R r^K ^r^A o/^R Z? /^A 



A^*,ALfA 



(4.20a) 
(4.20b) 
(4.20c) 



Let us now discuss two possible choices of F*, 



2. Out-scattering rate cutoff 



The simplest possibility is to choose 
F,,A = jFo,A ^ 0, 



(4.21) 



so that the cutoff-dependent distribution function defined 
via Eq. (|4.18p satisfies 



DZFoA ~ Fo.aZD + 2iAFo.A = 0. 



(4.22) 



In this case the QQ-componcnt of the inverse free prop- 
agator is chosen to be the following cutoff-dependent in- 
finitesimal regularization, 



(4.23) 



The term 2iAFo,A hi Eq. (|4.22p amounts to the following 
substitution for the time-derivative in the equations of 
motion for the distribution function, 



dt^dt + 2A. 



(4.24) 



The time-diagonal element of the non-interacting distri- 
bution function is then modified as follows, 



Fo,A(fc,i) =e-^^*Fo(fc,i) 



(4.25) 



whereas the cutoff-dependent non-interacting Keldysh 
Green function is now given by 



G^^A{k,t,t') = e-^^'+'">G^{k,t,t'). 



(4.26) 



The occupation numbers therefore decrease exponen- 
tially with rate A for large time, which justifies the name 
"out-scattering cutoff scheme." Because for A — > oo all 
propagators vanish in this scheme, the ERG fiow equa- 
tions should be integrated with the initial condition 

l™ ri"^,...o„ =0 ifn^4, (4.27) 

A—)- DO 

and the limit of T^a \ ^ ^ ^ is given by the bare inter ac- 
tiouEq. (2.51). 



3. Hybridization cutoff 

Alternatively, we may choose F^^a ~ Fq^a, so that the 
distribution function satisfies 



DZFo,A - Fo^aZD = 



(4.28) 



This equation agrees exactly with the cutoff-independent 
non-interacting kinetic equation (|2.44p . so that we may 
identify Fq.a = Fq with the cutoff-independent non- 
interacting distribution function. Obviously, this cutoff 
choice amounts to replacing the infinitesimal rj appearing 
in Eq. (|2.45p by the running cutoff A, so that 



(4.29) 



where Fq is the distribution function for infinitesimal 77, 
which is determined by the same equation as for A = 0. 
With this cutoff choice all propagators at non-equal times 
vanish for A — > 00. Explicitly, we obtain for the non- 
interacting 2x2 Green functions in flavor space in this 
cutoff scheme, 



G^^Aik,t,t')Fo{k,t')Z 
ZFo(fe,i)G^A(fc,i,0- 



(4.30) 



Because for large A —> 00 all propagators at non-equal 
times arc suppressed, each time integration in loops 
yields a factor of 1/A. For A — )■ 00 only the Hartree- 
Fock contribution to the self-energy survives because it 
depends only on the equal-time component of the Green 
function. The FRG flow equations in this cutoff scheme 
should therefore be integrated with the boundary condi- 
tion that for A — >■ 00 the irreducible self-energy is given 
by the self-consistent Hartree-Fock approximation. 

The finite value of the QQ-block of the inverse prop- 
agator can be considered to be a part of the Keldysh 
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self-energy, which in turn is related to an in-scattering 
ratei^ Compared to the out-scattering rate cutoff intro- 
duced before, this cutoff scheme contains both in and 
out scattering contributions, such that the bare distribu- 
tion function is cutoff independent. Essentially the same 
cutoff scheme has recently been proposed and tested by 
Jakobs, Pletyukhov and Schoelleri^i^^i^l In particular 
in Ref. [2^ they used this scheme to study stationary 
non-equilibrium states of the Anderson impurity model. 
In this case they identified the cutoff parameter A with 
the hybridization energy arising from the coupling to the 
bath of free electrons. Following their suggestion we shall 
therefore refer to this scheme as the "hybridization cutoff 
scheme." 



Substituting this approximation into Eq. (j4.20cp , we ob- 
tain for the Keldysh block of the single-scale propagator 
at equal times 

Gfik,t,t)^2i f diiGf(fc,t,ti)[FA(fc,ti)-F,,A(fc)] 

Jto 

xGiikM.t). (4.34) 

Obviously, this expression vanishes if the flowing distri- 
bution matrix Fp^ik^ti) approaches the equilibrium dis- 
tribution F*_A(fc)- 

C. Combining FRG flow equations with quantum 
kinetic equations 



4- Alternative cutoff schemes 

At this point it is not clear which cutoff choice is su- 
perior. By construction, both schemes do no violate 
causality for any value of the running cutoff A. More- 
over, to describe systems close to thermal equilibrium 
it might be important to require that in thermal equi- 
librium the fluctuation-dissipation theorem relating the 
Keldysh Green function to its retarded and advanced 
counter-parts is satisfied for any value of the running 
cutoff. Using the spectral representation of the Green 
functions, it is easy to show that for our model the 
fluctuation-dissipation theorem can be written as the fol- 
lowing relation between the Fourier transforms of the 2x2 
matrix Green functions in flavor space. 



G^(fc,cj) = [G^<(fc,w)-G^(fc,Lj)] 



1 



(4.31) 

While this relation is manifestly violated in the out- 
scattering cutoff scheme, in the hybridization cutoff 
scheme it remains valid in the non-interacting limit. 

Finally, let us point out that in certain situations other 
choices of the distribution matrix i^,.A in the QQ-block 
of the regularized inverse propagator given in Eq. (j4.15p 
may be advantageous. For example, for describing the 
approach to thermal equilibrium, it might be useful to 
identify i^,.A with the true equilibrium distribution. To 
see this, let us approximate the Keldysh Green function 
appearing in the Keldysh-block of the single-scale 
propagator (|4.20cP by the generalized Kadanoff-Baym 
ansatz (j3.24l) . which assumes that the non-interacting 
relation (j2.72l) remains approximately valid for the in- 
teracting system. Introducing the flowing distribution 
function in analogy with Eq. (|2.70|) , 



Ff,{k,t) = iZG'}^{k,t,t)Z, 



(4.32) 



the generalized Kadanoff-Baym ansatz p.24p can also be 
written as 



Gf{k,t,t')^Gf{k,t,t')Ff,{k,t')Z 
- ZFA{k,t)Gi{k,t,t'). 



The FRG flow equation (14.81) relates the derivative of 

the self-energy [EA]aia2 = ^A aia2 ^^^^ respect to the 
flow parameter A to the flowing Green function Ga and 
to the flowing effective interaction fI'^I, „ . Our final 
goal is to obtain a closed equation for the Keldysh block 
Ga of the Green function matrix at equal times (or al- 
ternatively the distribution function -Fa), from which we 
can extract the time evolution of the diagonal and off- 
diagonal distribution functions given in Eq. (|2.1al2.1bp . 
Therefore we have to solve the FRG flow equation (|4.8p 
simultaneously with the cutoff-dependent quantum ki- 
netic equation, which can be derived analogously to 
Sec, mil from the cutoff-dependent Dyson equation. 



(4.35) 



The cutoff-dependent kinetic equation for the Keldysh 
block can be derived in the same way as in Sec. IIIIl and 
we thus obtain for the equal-time Keldysh Green func- 
tion. 



idtG^ik, t, t) - MA,kG^ (fc, t, t) ~ G^ (fe, t, t)Ml^ 



K I 



f dh[zj:f{k,t,ti)Gi{k,h,t) 

-G^(fc,t,ii)Sf(fc,ti,t)Z] 

/ dh[zj:'l{k,t,h)G'l{k,h,t) 

Jtn 



-Gf{k,t,ti)^i{k,ti,t)Z],{4.M) 



(4.33) 



where Ma^u is a cutoff dependent deformation of the ma- 
trix Mfc defined in Eq. (|2.42p . The explicit form of Mam 
depends on the cutoff scheme. For the out-scattering 
cutoff scheme, it follows from Eq. (|4.22p and (|4.32p that 
Mam = AJ^k — *A/. For the hybridization cutoff scheme, 
Eqs. and imply that Mam = ^^k, which 

is identical to Eq. (|2.42p . The general form of the ki- 
netic equation (|4.36p is of course similar to Eq. p.23p . 
except that now all Green functions and self-energies de- 
pend on the cutoff parameter A. Together with the FRG 
flow equation (j4.8l) , this equation forms a system of cou- 
pled first-order partial intcgro-diffcrential equations with 
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two independent variables t and A, which have to be 
solved simultaneously. Because the flow equation (|4.8p 
depends on the effective interaction which satisfies the 
flow equation (|4.9p . the simplest truncation is to neglect 
the flow of the interaction. However, the resulting sys- 
tem of kinetic and flow equations given by Eqs. (|4.36p and 
(|4.8p is not closed because the flow and the kinetic equa- 
tion contains integrals involving the two-time Keldysh 
Green function. To reduce the complexity and to close 
the system of equations, the usual approximation strate- 
gies of quantum kinetics can now be made. For exam- 
ple, on the right-hand side of the quantum kinetic equa- 
tion (|4.36p . one could express the Keldysh Green func- 
tion for non-equal times in terms of the corresponding 
equal-time Keldysh Green function using the generalized 
Kadanoff-Baym ansatz (|3.24p . Further simplifying ap- 
proximations such as the Markov-approximation where 
the time-arguments of all Keldysh Green functions on 
the right-hand side of Eq. (|4.36p are replaced by the ex- 
ternal time t might also be useful. For the solution to be 
unique, we have to further specify the boundary condi- 
tions. For our system of kinetic and flow equation it is 
sufficient to define the distribution i^A(io) at the initial 
time to for arbitrary cutoff A, and the self-energy Sao (*) 
at the initial cutoff scale Aq for arbitrary time t. We shall 
illustrate this choice of boundary conditions and the ap- 
proximations mentioned above in the following section 
within the framework of a simple exactly solvable toy 
model. 



For notational simplicity, we redefine again^" e — ^ — 
e. This simplified model describes a single anhar- 
monic quantum mechanical oscillator subject to a time- 
dependent external field which creates and annihilates 
pairs of excitations. Although this toy model does not 
describe relaxation and dissipation processes, it does cap- 
tures some aspects of the physics of parametric resonance 
in dipolar ferromagnetsi^ 

The non-equilibrium dynamics of the Hamiltonian 
(|5.2p can be easily determined numerically by directly 
solving the time-dependent Schrodinger equation. Ex- 
panding the time-dependent states \ip{t)) of the Hilbert 
space associated with Eq. (|5.2p in the basis of eigenstates 
\n) of the particle number operator a^a, 



(5.3) 



n=Q 



the time-dependent Schrodinger equation assumes the 
form 

r u 1 
idti>n{t) = en + -n(n - 1) Tpnit) 

, M 

2 



y/n{n - l)V',i-2(i) 
+ ^in + 2)in+l)i'^+2it)\. (5.4) 



This system of equations is easily solved numerically. 
From the solution we may construct the normal and 
anomalous distribution functions, 



V. EXACTLY SOLVABLE TOY MODEL 

Although the functional renormalization group ap- 
proach for bosons out of equilibrium developed in Sec. lIVI 
is rather general, at this point it is perhaps not so clear 
whether this approach is useful in practice to calculate 
the non-equilibrium time evolution of interacting bosons. 
One obvious problem which we have not addressed so far 
is that truncation strategies of the formally exact hierar- 
chy of ERG flow equations have to be constructed which 
correctly describe the long-time asymptotics. 

As a first step in this direction, we shall consider in 
this section a simplified version of our boson Hamiltonian 
p.ip which is obtained by retaining only the operators 
associated with the fc = mode. Setting ak=o = a, 
efc=o = 7fc=o = 7 and C/(0, 0; 0, 0)/F = u, our boson 
Hamiltonian (jl.ip thus reduces to the following bosonic 
"toy model" Hamiltonian, 



nit) = ea^a + ^ ['je-"^''*a'<a'< + 7*e*'^°* 



aa\ 



" t t 
-|- —a'a'aa. 



(5.1) 



In the rotating reference frame Eq. (|5.ip becomes 



n 



LOq 



^ ^ a'^a + ^ [a^'a^ + aa] 



-a^a^aa. (5.2) 



nit) = (mWam)) = E "IV'n Wl"' (s-s) 

ri=0 

p{t) = (^(i)|aa|^(i)) 

OO 

= E \/(n + 2)(7i + l)v:(t)Vn+2(i). (5.6) 

n=0 

We have prepared the coefficients at the initial time as 

'ipnito) = Sn,l, (5.7) 

SO that the normal and anomalous pair-correlators have 
the initial values n(0) = 1 and p{0) — 0. With this 
choice all other correlators vanish at the initial time. We 
have solved the Schrodinger equation (|5.4p numerically 
by both integrating it directly and by calculating the ma- 
trix exponential exp[— i'H(t — to)] and using 



E 

3=0 



^3 (to), 



(5.8) 



where the Hamiltonian Ji has the matrix elements 
1) 



^ nrn 



u 

-n{n 



y/in + 2){n+l)Sn+2,. 



l)'5?i-2,m 

(5.9) 



We found identical results with both methods. A total 
number of nmax — 20 basis coefficients was sufficient for 
convergence. 
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A. Time-dependent Hartree-Fock approximation 

As a reference, let is briefly discuss the self-consistent 
Hartree-Fock approximation for our toy model. In the 
context of parametric resonance of magnons in yttrium- 
iron garnet, this approximation is also referred to as "S- 
thcory."— i2£ Within this approximation the self-energy 
matrix is time-diagonal. 



(5.10) 



With the help of the symmetrized interaction vertex de- 
fined in Eq. (j2.5ip we may write 

o'icr2 A1A2 

= ^ E E u^'J:'''G',i'.i it, t)- (5-11) 

RecaU that by definition E^c- ^ j^fl j^CQ ^ ^A^ 
so that we obtain for the time-diagonal elements of the 
retarded and advanced self-energy, 



p*it) 2n{t) + l 
2n(t) + 1 p{t) 



(5.12) 



The Keldysh component of the self-energy vanishes in 
this approximation. 



(5.13) 



Actually, there is an additional time- independent in- 
teraction correction —u to the normal component of 
the advanced and retarded self-energy which arises from 
the symmetrization of the Hamiltonian, as discussed in 
Sec. Ill CI According to Eq. (|2.48p this contribution 
simply leads to a constant shift —u in the energy in 
Eq. (|2.42p . Taking this shift into account, we find that 
our kinetic equation (|3.23p reduces to the following 2x2 
matrix equation. 



idtF{t) ^ -M^{t)F{t) - F{t)M{t), 



where 



M{t) = M + ZY.i{t) = 



with 



and 



M = 



e{t) 7(0 



l7l 



H7I -(£-") 



e(t) = e + 2un{t), 
lit) = h\+up{t). 



(5.14) 



(5.15) 



(5.16) 



(5.17a) 
(5.17b) 



Recall that according to Eq. p.70p the 2x2 distribution 
matrix is given by 



F(.).^ZG-(M)^=t2n(r);i 



-2p*{t) 2n{t) + l 
-2p{t) 



(5.18) 



At this level of approximation the kinetic equation (|5.14p 
has the same structure as the corresponding equation 
(|2.7ip in the absence of interactions. From Eqs. (|5.18p 
and (|5.14p we obtain the following kinetic equations for 
the diagonal and off-diagonal distribution functions^ 

idtn{t) = -i*{t)p{t) + l{t)p*{t), (5.19a) 
idtp{t) = 2e{t)p{t) + i{t)[2n{t) + 1]. (5.19b) 

In Figs. [T] and [2] we compare the numerical solution 
of these equations with the exact result obtained from 
Eqs. (j5.4H5.6p . and with the time evolution in the non- 
interacting limit. Because our simple toy model does 
not account for damping and dissipative eflFects, the time 
evolution is purely oscillatory. However, the true oscil- 
lation period lies between the non-interacting oscillation 
period Tq = -k j p. = -k j \/ — |7p and the smaller oscilla- 
tion period predicted by the self-consistent Hartree-Fock 
approximation. A similar phenomenon is also observed 
for the oscillation amplitudes. As expected the devia- 
tion between the three curves increases with increasing 
interaction strength. We thus conclude that the Hartree- 
Fock approximation is only capable for moderate interac- 
tion strength and short times as illustrated in the middle 
panel of Figs. [l]and[2l In this case the Hartree-Fock re- 
sult up to times of the order To/4 is reliable. However, in 
this regime the free time evolution is also fairly accurate. 



B. Kinetic equation with self-energy up to second 
order 



Let us consider again the quantum kinetic equation for 
the Keldysh Green function G^{t,t'), which for our toy 
model can be obtained by simply omitting the momen- 
tum label in Eqs. p. lip and (|3.12p . Substituting on the 
right-hand side of this equation the self-energies up to 
second order in the interaction given in Eq. (|5.1ip (first 
order self-energy) and in appendix C (second order self- 
energy), we obtain an equation of motion for the two- 
time Keldysh Green function G^{t,t'). Together with 
the corresponding retarded and advanced components of 
the Dyson equation, this equation forms a closed system 
of partial differential equations, which can in principle be 
solved numerically. To simplify the numerics, we will fo- 
cus here only on the evolution equation for the equal-time 
Keldysh Green function G^{t,t) which can be obtained 
from the kinetic equation p.23p by omitting the momen- 
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FIG. 1. (Color online) Time evolution of the diagonal distri- 
bution function for initial conditions n(0) = 1 and p(0) = 0. 
We have chosen |7|/e = 0.1 and three different values of the 
interaction strength: u/e = 0.025 (top), u/e = 0.1 (middle), 
and u/e = 0.25 (bottom). The frequency = — |7p de- 
termines the oscillation period in the non-interacting limit. 
We compare the result of the self-consistent Hartree-Fock ap- 
proximation (solid line) with the exact solution (dashed line) 
and the time evolution in the non-interacting limit (dotted 
line) . 



turn labels, 

idtG^{t,t) - MG^(t,t) - G^(t,t)M^ 



f dti[ZS^(i,ti)G^(ii,i) 

J to 

-G^(i,ti)I]^(<i,t)Z] 
[ dti[ZS^(t,ti)G^'^(ti,t) 

J to 

-G^(i,ti)I]^(ii,t)Z] 
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FIG. 2. (Color online) Time evolution of the off-diagonal 
distribution function for increasing interaction strength (from 
top to bottom). The parameters and initial conditions are 
the same as in Fig. [T] We compare the result of the self- 
consistent Hartree-Fock approximation (solid line) with the 
exact solution (dashed line) and the time evolution in the 
non-interacting limit (dotted line). 



Since the two-time function G^{t,t') appears again on 
the right hand side of this equation let us make three 
additional standard approximations to close the system 
of equations: 

1. Generalized Kadanoff-Baym ansatz: As discussed 
in Sec. IIII Bl with the help of the generalized 
KadanofF-Baym ansatz (|3.24p we may derive a 
closed integral equation for the equal-time Keldysh 
Green function G^{t,t). For our toy model, the 
generalized KadanofF-Baym ansatz reads 



(5.20) 



G^{t,t') « ~i[G^{t,t')ZG^{t',t') 
-G^it,t)ZG^it,t')]. 



(5.21) 
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This ansatz amounts to approximating the distri- 
bution matrix in the colHsion integrals by its diag- 
onal elements, 



F{t,t') 



5{t 
5{t 



t')F{t) 

t')iZG^{t,t)Z. 



(5.22) 



2. Markov approximation: To reduce the integro- 
differential equation for the equal-time Keldysh 
Green function to an ordinary differential equation, 
we replace under the integral 



diiG'"(ti,ti)... G'"(i,0 / dti.... (5.23) 



3. Free advanced and retarded propagators: Finally, 
we neglect the self-energy corrections of the ad- 
vanced and retarded propagators in Eq. ()5.20|) and 
thus replace G^{t, t') and G^{t, t') by the free prop- 
agators, which can be obtained by omitting the mo- 
mentum labels in Eqs. (|2.55l2.60p . 

After these approximations, the collision integrals in 
Eq. (|5.20|) can be calculated analytically and the non- 
equilibrium distribution functions are easily obtained by 
numerically solving a system of two coupled ordinary dif- 
ferential equations. For the numerical solution, the time 
grid was chosen equally spaced with Ate = 1.3 x 10^^ 
and the differential equations were solved using a fourth- 
order Runge-Kutta algorithmi^ 

The result for the same parameters and initial con- 
ditions as in Figs. [1] and [2] is shown in Figs. |3] and 
m Obviously, for moderate interaction strength the in- 
clusion of second-order corrections indeed improves the 
agreement with the exact solution up to times of order 
To = t: I n = 7r/-\/e2 — |7p. However, for stronger in- 
teractions and for times exceeding Tq, the solution of 
the kinetic equation with second-order corrections to the 
self-energy disagrees even more drastically from the exact 
solution than the time-dependent Hartree-Fock approx- 
imation shown in Figs. [1] and [5J In addition we found 
secular behavior and unphysical divergences of the pair 
correlators for long times (not shown in Figs. [3] and |4]). 
By numerically solving the kinetic equation without mak- 
ing the above approximations (see appendix D), we have 
checked that the strong disagreement of the time evolu- 
tion beyond one oscillation period Tq with the exact re- 
sult is not an artifact of the Kadanoff-Baym ansatz, the 
Markov approximation or the neglected renormalization 
of the retarded and advanced propagators. 



C. First order truncation of the FRG hierarchy 
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FIG. 3. (Color online) Time evolution of the diagonal distri- 
bution function of the toy model. The parameters and initial 
conditions are the same as in Fig. [1] The solid line is the so- 
lution of the kinetic equation with second order self-energies, 
simplified using the generalized Kadanoff-Baym ansatz and 
the Markov approximation. For comparison we also show the 
exact solution (dashed line). 



much better results for the time evolution than the pre- 
vious two approximations. A similar truncation has also 
been made by Gezzi et ali^ in their FRG study of sta- 
tionary non-equilibrium states of the Anderson impurity 
model. In the exact FRG flow equation (|4.8p for the 
self-energy we replace the flowing effective interaction by 
the bare interaction (recall that for our toy model the 
collective labels ai represents (t^, Ai,CTi)), 



Wc now show that a very simple truncation of the non- 
equilibrium FRG flow equation for the self-energy where 
the flow of the effective interaction is neglected leads to 



" A,QiQ;2Ct3Ct4 



6{ti - t2)5{t2 - t3)5{t3 - ti)U^l 



Ai A2 A3 A4 

(5.24) 
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FIG. 4. (Color online) Time evolution of the off-diagonal dis- 
tribution function for increasing interaction strength (from 
top to bottom) of the toy model. The parameters and initial 
conditions are the same as in Fig. [1] The solid line is the so- 
lution of the kinetic equation with second order self-energies, 
simplified using the generalized Kadanoff-Baym ansatz and 
the Markov approximation. The dashed line is again the ex- 
act solution. 



where up to permutation of the indices the symmetrized 
bare interaction is given by (see Eq. (|2.5ip ) 

jjCQCC ^ jjCQQQ ^ jjCCCQ ^ jjQQCQ ^ 

In this approximation the two-point function is time- 
diagonal, 

^fX't'=5{t-t')Y.i%.it), (5.26) 
where the self-energies satisfy the FRG flow equation 

9AEi^:,,(t) = ^ E E u5ll^^lGi%^_{t,t). (5.27) 

AlfTl A2Cr2 



With the bare interaction given by Eq. (j5.25p . this leads 
to the FRG flow equations for the retarded (AA' = QC) 
and advanced (AA' = CQ) self-energies, 



aaii.t) 



(5.28) 



For the Keldysh component of the self-energy corre- 
sponding to AA' — QQ we obtain 



it,t) 



(5.29a) 
(5.29b) 



From the definitions (|4.20al4.20bp of the retarded and 
advanced components of the single-scale propagators it 
is easy to see that at equal times G^(^' t) = = G^{t, t), 
so that within our truncation the right-hand sides of the 
flow-equations (|5.29al5.29bp for the Keldysh self-energy 
vanish. Because the initial Keldysh self-energy is zero, it 
remains zero during the entire RG flow within our trun- 
cation. 

At this point we specify our cutoff procedure. It 
turns out that from the two cutoff schemes discussed in 
Sec. IIVBI the out-scattering rate scheme described in 
Sec. IIVB"2] is superior. Recall that this scheme amounts 
to setting F^^a = in Eq. (j4.20cp . For our toy model 
the Keldysh component of the single-scale propagator at 
equal times is then given by the following 2x2 matrix 
equation. 



Gf{t,t) =i dti 

Jto 



Gf{tM)ZGf{h,t) 
G^{t,h)ZGi{h,t) 



(5.30) 



Note that this equation still contains memory effects. We 
simplify Eq. (|5.30|) using the same approximations as in 
the previous subsection: First of all, we use the general- 
ized Kadanoff-Baym ansatz (|5.2ip to express the Keldysh 
Green functions on the right-hand side of Eq. (|5.30p in 
terms of the corresponding equal-time Green function. 
Introducing the cutoff-dependent distribution function 



FA{t)=tZGf{t,t)Z, 



(5.31) 



we obtain 



G^{t,t)^2t dhG^{t,h)FA{h)Gi{h,t). (5.32) 



Next, we replace FA(ii) — >• FA{t) under the integral sign 
(Markov approximation) . Substituting the advanced and 
retarded propagators by their free counterparts (which 
can be obtained from Eq. (|2.55l2.60p by omitting the mo- 
mentum labels) we finally arrive at the following simple 
expression for the Keldysh component of the single-scale 
propagator in the out-scattering rate cutoff scheme, 

G^it,t)^2i f dhG^j,{t,h)FA{t)G^^A{tut). (5.33) 
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At this point we have arrived at a system of two cou- 
pled partial differential equations (PDE) for the cutoff- 
dependent distribution matrix F^it) and the self-energy 
matrix T,\(t). The former contains the normal {n/^{t)) 
and anomalous {pA{t)) distributions as in Eq. (|5.18p . 

The time evolution of this distribution matrix is deter- 
mined by a kinetic equation which is formally identical 
to the corresponding kinetic equation (|5.14p within time- 
dependent Hartrcc-Fock approximation, 

idtFAit) = -Ml{t)FA{t) - FA{t)MA{t). (5.35) 

The cutoff-dependent matrix M^it) is 

MA{t) = M - iAI + ZT.A{t) 

^ ( e-u-iK + Y.A^aa{t) |7l + SA.aa(<) A 

^ -|7| -SA,aa(i) ~[t-U + iK + ^A,a-a{i)] ) ' 

(5.36) 

and depends on the bare matrix M defined in Eq. (j5.16p , 
on the cutoff parameter A, and on the cutoff-dependent 
self-energy SA(i)- The flowing self-energy matrix satisfies 

where the matrix G^{t,t) is given by Eq. (|5.33p . 

Mathematically, the problem is now reduced to the so- 
lution of a system of first order partial differential equa- 
tions in two independent variables t and A. To illustrate 
the structure more clearly, we rewrite the system ()5.35|) 
and (|5.37p as 

dtFAit) - A{FA{t),J:Ait),A), (5.38a) 
dA^Ait)^B{FA{t),A,t), (5.38b) 

where the explicit form of the matrix functions A and 
B follows from the right-hand sides of Eq. (|5.35p and 
(|5.37p . Note, that the system is not fully symmetric in 
the variables A and t, because the flow equation con- 
tains causal memory integrals over the time t. Without 
Markov approximation, the right-hand side of the flow 
equation (and in higher-order truncations also the ki- 
netic equation) is a functional of the distribution matrix 
FA{t) and depends on the distribution matrix at earlier 
times. Using the Markov approximation the distribution 
matrix can be pulled out of the integral and the func- 
tional B reduces to an ordinary function of the distribu- 
tion FA{t). To deflne the solution of Eqs. (|5.38al5.38bp 
uniquely we note that the boundary conditions fix the 
distribution matrix i^A(to) the initial time to and ar- 
bitrary cutoff A, and the self-energy matrix EAo(i) at 
the initial cutoff Aq and arbitrary time t. In fact, within 
our truncation, the boundary condition for the distribu- 
tion matrix is Ja(^o) = F{to). Since for a large cutoff 




FIG. 5. (Color online) Illustration of our approach to solve 
the system (|5.38a|) and ()5.38b|) of partial differential equa- 
tions. The kinetic equation (|5.38ap describes the propagation 
of FA{t) in t, where the flow equation (|5.38b[) describes the 
propagation of SA(t) in A-direction. The boundary conditions 
define the distribution function FA{to) at the initial time to 
(dashed line) and the self-energy Eaq (t) at the initial Aq (dot- 
ted line). 

all one-particle irreducible vertices (with the exception 
of F^^)) vanishes due to Eq. (|4.27p . the boundary condi- 
tion for the self-energy matrix at sufficiently large initial 
cutoff Aq is SAo(i) = 0. The standard method of deal- 
ing with this kind of first order PDEs is the method of 
characteristics!^ However, in our case the characteristic 
curves coincide with the curves where the boundary con- 
ditions are specified, so that the standard procedure is 
not applicable. Nevertheless, it is easy to see that the 
solution with the proper boundary conditions can be ob- 
tained by means of the following algorithm: We first note 
that the kinetic equation (|5.38ap describes the propaga- 
tion of FA{t) in t, and that the flow equation (j5.38bp 
gives the propagation of SA(t) in A direction, as illus- 
trated in Fig. m Solving the kinetic equation (|5.38ap for 
an inflnitesimally small time step dt, the resulting dis- 
tribution function at t + dt can be used to integrate the 
flow equation (j5.38bp at flxcd t + dt over A. Repeating 
these two steps allows to obtain the solution of Fa (t) and 
SA(i) in the entire (i, A)-planc. 

In the following we explain our approach to numeri- 
cally solve the coupled set of flrst order partial differ- 
ential equations. We focus on the out-scattering cutoff 
scheme but generalizations to other cases are straightfor- 
ward. Wc consider a discretization of the two variables t 
and A in the form 

t^t,ne{to,...,tM~i}, (5.39a) 
A^ A„ G {Ao,...,AAr_i}, (5.39b) 

with m G {0, . . . , A/ - 1} and n G {0, . . . , iV - 1}. The 
discretized grid points are ordered as i„ < tn+i and 
An > A„+i. Both grids do not need to be equally spaced. 
Moreover, the number of points M and A'^ can be chosen 
arbitrarily. The discretized functions arc written as 

FA„{tm) ^ Fnm , SA„(tm) = S„m, (5.40) 
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and 



n m ■ 



(5.41) 



The derivatives are approximated by first-order finite- 
difference expressions, 



-^n m+1 


— F 


^?7l + l 


— f ' 






An+1 


- A„ 



(5.42a) 
(5.42b) 



The discretized version of tlie kinetic equation (|5.14p tlien 
follows as 

(5.43) 

with the time-dependent coefficient matrix 



Mnm =M - iAnI + ZY.nm- (5.44) 



In the same way the discretized flow equation ()5.37|) for 
the self-energy is 

+ IU[K+, - A„) ( ^^^^ra.ra-a G^^,a-a \ (5 45) 
\ nm,aa 2 nm^aa / 

According to Eq. (|5.33p the single-scale propagator 
on the right-hand side is defined as 



to 



(5.46) 

From the structure of the discretized equations it is obvi- 
ous that causality is preserved since each time step can be 
calculated from the previous ones and does not depend 
on quantities at later times. Starting from the initial 
values which specify the distribution matrix i^„o with 
n g {0, . . . , A'^ — 1} and the self-energy matrix Som with 
m S {0, . . . , M — 1} on the boundaries, one can obtain 
the entire solution by stepwise propagating in t- and in 
A-direction in terms of basic Euler steps. One Euler step 
from tm to tm+i contains two parts: First, with the so- 
lution Fnm where n G {0, . . . , — 1} from the previous 
step and the initial self-energy Eom on the boundary, 
the flow equation (|5.45p at fixed time can be inte- 
grated in sub-steps from Aq to A^r-i to obtain 
on all points n G {0, . . . , — 1}. Next, by using the ki- 
netic equation (|5.43p . Fnm+i can be derived from Fnm 
and Snm- This completes one basic Euler step since the 
distribution function _F„,„+i is now known at time i,„+i. 
Repeating this procedure M-times yields the full solution 
up to time <A/-i- Numerically, the first-order finite differ- 
ence derivatives are not accurate enough unless the grid 
spacing becomes very small which is not feasible in prac- 
tice. Therefore a fourth-order Runge-Kutta method for 
the propagation in t- and the second-order Heun method 



for the propagation in A-directioE^^ is used. One Runge- 
Kutta step from Fnm to Fnm+i consists of four Euler 
steps of the form described above. The integral (|5.46p 
was solved analytically using the free retarded and ad- 
vanced Green functions (given by p. 551 I2.60p without 
k dependence). The time grid was chosen similar to 
the Hartrec-Fock and the second-order case described in 
Sec. IV Bl The A-grid ranges between Ao/e = 8.1 and 
A299/e = 2.1 X 10"'' and was adjusted in such a way that 
for A/e < 1, the resolution of the grid spacing was in- 
creased to take into account the higher curvature of the 
self-energy in this region. 

The result for the FRG approach with the out- 
scattering cutoff scheme for the same parameters and ini- 
tial conditions as in perturbation theory (compare Figs. 
[1] [5] and [21 d]) is shown in Figs. [B] and [71 For the oscil- 
lation period of the pair correlators, the FRG treatment 
clearly improves the results compared to perturbative ap- 
proaches. Up to time To — TT / fi — TT / e'^ — |7p the pe- 
riod of the oscillation is nearly identical with the exact 
result. Even after longer times of the order 4To (mid- 
dle panel) the deviation from the exact solution is small. 
The oscillation is regular and we found no secular behav- 
ior even at long times. However, the amplitude of the 
pair correlators is underestimated and is comparable to 
the perturbative mean-field result shown in Figs. [1] and 
[21 In contrast, with the hybridization cutoff scheme, we 
were not able to obtain any reasonable results for the 
pair-correlator dynamics. This suggests that in practice 
the out-scattering cutoff scheme works better than the 
hybridization cutoff scheme. 



VI. SUMMARY AND CONCLUSIONS 

We have developed a real-time functional renormaliza- 
tion group (FRG) approach to calculate the time evolu- 
tion of interacting bosons out of equilibrium. To be spe- 
cific, we have developed our formalism in the context of 
the interacting time-dependent boson Hamiltonian (|l.ip 
describing the non-equilibrium dynamics of magnons in 
dipolar magnets such as yttrium-iron garnet^ subject 
to an oscillating microwave field i^^i^° To take into ac- 
count the off-diagonal correlations inherent in this model 
we have introduced an efffcient matrix notation which 
facilitates the derivation of quantum kinetic equations 
for both the normal and anomalous components of the 
Green functions in the Keldysh formalism. We have also 
extended the generalized Kadanoff-Baym ansata^ to in- 
clude both diagonal and off-diagonal correlations on the 
same footing. 

In our FRG approach the time evolution of the di- 
agonal and off-diagonal distribution functions is ob- 
tained by solving a quantum kinetic equation with cutoff- 
dependent collision integrals simultaneously with a renor- 
malization group flow equation for cutoff-dependent non- 
equilibrium self-energies appearing in the collision in- 
tegrals. To implement this procedure, we proposed a 
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FIG. 6. (Color online) The solid lines are our FRG results for 
the diagonal distribution function for increasing interaction 
strength (from top to bottom). The parameters and initial 
conditions are the same as in Figs. [T]and|3] For comparison 
we also show the exact solution (dashed line). 
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for the off-diagonal distribution function for increasing inter- 
action strength (from top to bottom). The parameters and 
initial conditions are the same as before. The dashed line is 
again the exact solution. 



new cutoff scheme where the infinitesimal imaginary part 
defining the boundary conditions of the inverse advanced 
and retarded propagators is replaced by a finite scale 
acting as a running cutoff. We have called this cutoff 
procedure out-scattering rate cutoff scheme because the 
cutoff-dependent imaginary parts in the retarded and 
advanced propagators lead to an exponential decay of 
the occupation numbers. In principle one can also re- 
place the infinitesimal imaginary part appearing in the 
Keldysh component of the inverse free propagator by a 
cutoff-dependent finite quantity which leads to the hy- 
bridization cutoff scheme proposed by Jakobs et al.^ 
For the toy model we presented evidence that it is better 
to work with the out-scattering cutoff scheme, keeping 
the Keldysh component of the inverse free propagator 



infinitesimal. 

We have explicitly tested our FRG approach for a sim- 
plified toy model which is obtained from the Hamilto- 
nian (|l.ip by retaining only a single momentum mode. 
Although this simplified model does not contain damp- 
ing and dissipative effects, it does describe some aspects 
of the magnon dynamics in yttrium-iron garneti^ Since 
the non-equilibrium time evolution of our toy model can 
be obtained exactly by direct numerical integration of 
the time-dependent Schrodinger equation, our toy model 
allows us to test the quality of various approximations. 
Specifically, we have studied the following approxima- 
tions: 

1. Self-consistent Hartrec-Fock approximation, which 
is also called S-thcory in the context of non- 
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equilibrium dynamics of magnonsi22i2£ 

2. A pcrturbative approach based on the calculation 
of the non-equilibrium self-energies to second order 
in the interaction, in combination with the gener- 
alized Kadanoff-Baym ansatz and the Markov ap- 
proximation. 

3. A FRG approach based on the simultaneous solu- 
tion of a coupled system of kinetic equations and 
renormalization group flow equations for the scale- 
dependent self-energies, using a simple truncation 
of the FRG flow equations for the non-equilibrium 
self-energies where the flow of the interaction is ne- 
glected. 

For each approach we have calculated the time depen- 
dence of the normal and anomalous distribution function 
for some representative values of the interaction and com- 
pared the result with the exact solution. It turns out that 
the first two approaches do not give reliable predictions 
for the time evolution beyond one oscillation period. Al- 
though inclusion of second order self-energy corrections 
somewhat improves the agreement for short times, the 
time dependence beyond a single oscillation period dis- 
agrees even more strongly with the exact solution than 
the prediction of the self-consistent Hartree-Fock approx- 
imation. In fact, in Appendix D we shall show that the 
numerical solution of the two-time quantum kinetic equa- 
tions with second order self-energies does not lead to a 
better agreement with the exact solution of the toy model 
than calculations which in addition rely on the Kadanoff- 
Baym ansatz and the Markov-approximation. The per- 
turbative approaches are therefore not able to reproduce 
the real-time dynamics of our toy model and do not al- 
low for systematic improvements. The failure of pertur- 
bation theory to predict the long-time behavior of cor- 
relation functions is not unexpectedJ^ In contrast, our 
simple truncation of the FRG flow equations in combina- 
tion with the out-scattering cutoff scheme leads to quite 
good agreement with the exact solution over many oscil- 
lation periods. Note, however, that our FRG approach 
is numerically more costly than the other two methods, 
because one has to solve a coupled system of partial dif- 
ferential equations in two independent variables, the time 
t and the cutoff-parameter A. Moreover, due to our sim- 
ple truncation of the FRG flow equations the oscillation 
amplitudes are still underestimated. Nevertheless, from 
all approximation strategies we have tested our first order 
FRG approach with out-scattering cutoff scheme clearly 
gives the most satisfactory results for the time-evolution 
of our toy model. One should keep in mind, however, 
that the toy model does not have any intrinsic dissipa- 
tion and therefore does not relax towards a stationary 
state at long times. In fact, we expect that for mod- 
els with intrinsic dissipation other cutoff schemes such as 
the hybridization cutofS^i^l discussed in Sec. lIVBj is su- 
perior, because the hybridization cutoff retains the bal- 
ance between in-scattering and out-scattering terms in 



the collision integral which is crucial to describe the re- 
laxation towards a stationary state. The hybridization 
cutoff scheme also preserves the fluctuation-dissipation 
theorem during the entire flow2^ which is advantageous 
close to thermal equilibrium. 

Our work can be extended in several directions: First 
of all, it should be interesting to use our FRG approach 
to calculate the time evolution of infinite or open quan- 
tum systems which exhibit relaxation and dissipative pro- 
cesses. We expect that for such systems standard ap- 
proximations such as the Kadanoff-Baym ansatz or the 
Markov approximation may have different regimes of va- 
lidity than for our toy model. Recall that for weakly 
correlated systems such as semiconductors the Kadanoff 
Baym ansatz has been shown to be quite useful and 
accurate!^ It should also be interesting to extend our 
study of the toy model to the regime of strong pump- 
ing where the original vacuum state is unstablci^ More- 
over, it would be even more interesting to apply our non- 
perturbative FRG method to study the non-equilibrium 
dynamics of the time-dependent boson Hamiltonian (jl.ip 
in the regime of strong pumping. It is well know n^^i'^° 
that for sufficiently large values of the pumping parame- 
ter 7fc the system exhibits the phenomenon of paramet- 
ric resonance. The magnon operators acquire than finite 
expectation values and the system approaches a non- 
trivial time-independent non-equilibrium state which is 
dominated by interactions! • Although this state 
has been studied at the level of time-dependent Hartree- 
Fock approximationSSi^S (S-theory) it would be interest- 
ing to describe the time evolution into stationary non- 
equilibrium states non-perturbatively, and check if the 
states exhibit non-thermal scaling properties as predicted 
in Ref. [l^. 



ACKNOWLEDGMENTS 

The authors thank J. Berges and J. Hick for useful dis- 
cussions. Financial support by SFB/TRR49, FOR 723, 
the CNRS, and the Humboldt foundation is gratefully 
acknowledged. 



APPENDIX A: TRANSFORMATION TO THE 
ROTATING REFERENCE FRAME 

In order to simplify the calculations, it is useful to re- 
move the explicit time-dependence from the Hamiltonian 
7i{t) in Eq. (jl.ip . This can be achieved by means of a 
unitary transformation to the rotating reference frame, 
as discussed in Sec. HI In the rotating reference frame 
the Hamiltonian does not explicitly depend on time, see 
Eqs. (|1.6p . To distinguish quantities in the original- and 
the corresponding rotating frame, we put in this ap- 
pendix an extra tilde over Green functions in the rotating 
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frame. Introducing the unitary 2x2 matrix 



Uk{t) 
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(wot-ipfc) 



(Al) 



the relations between the elements of the matrix Green- 
functions G", G^, G^ defined in Eqs. (12.8aH2.8c|) and 
the corresponding quantities in the rotating reference 
frame (with the tilde) is 

G^(fe,t,t') = Uk{t)G''{k,t,t')Uu{t'), (A2) 

where X = R,A, K, and we have defined the 2x2- 
matriccs in fiavor space, 



G''{k,t,t') = [G'']kt,-kt'. 
Introducing the diagonal matrix 



(A3) 



(A4) 



with matrix elements [U]kt,k't' = Sk.k'S{t — t')Uk(t), we 
can rewrite Eq. (jA2[) in the more compact form 



G' 



UG^U 



(A5) 



The distribution function matrix F defined via Eq. (|2.15p 

is related to its counterpart F in the rotating reference 
frame via 



F ^ U^'FU''. 



(A6) 



Taking matrix elements in the time-labels and using the 
fact that in the non-interacting limit the distribution 
function matrix is time-diagonal (see Eqs. ()2.68p ). the 
relation (|A6p implies the 2x2 matrix equation, 



(A7) 



Foik,t)^ul{t)Fo{k,t)ulit), 
where we have used the fact that 

Ulit) = Z^Uk{t)Z = ZUkit)Z^ 



(AS) 



so that by definition 

Gii,^Gf,^ + G^/. (B2) 

Note that the above functions are 2x2 matrices in flavor 
space. Acting with (G^)"^ from the left on Eq. (|Bla[) 
and using the left Dyson equation (|3.2a|) for the retarded 
Green function we obtain 



[iG'')-'G''%t' ^ ~id{t-t')ZG^t, 



+ e{t-t') [[{G§)-'G%t' - dhi:lGf^,, 

(B3) 

A similar relation can be derived for the advanced com- 
ponent of the Keldysh Green function, 

[G''^iG^)-%, ^ iS{t ^ t')Gft,Z 

+ e(t' - 1) (^[G'HG^r'W J' dt,Gf,X.t^ ■ 

(B4) 

Using the Keldysh components of the left and right 
Dyson equations given in Eqs. p.2cp and (|3.4cp . the ac- 
tion of the free inverse propagators on G^ can be written 
as 

[{G§r'G%t' = f dhJ:lGf^,, + f dhJ:f,^Gt„ 



(B5) 



(B6) 



Substituting these expressions into Eqs. (|B3p and (|B4p 
and solving for G^^ and we obtain 

[G''%t. =-iG^t,ZG^t'+Q{t-t') 

X fdt.f dt2G^JJl]^,,G^;,+S,^,^Git,,], 

Jt' J -oo 

(B7) 



APPENDIX B: GENERALIZED 
KADANOFF-BAYM ANSATZ 

The generalized Kadanoff-Baym ansatz (GKBA) is an 
approximate relation between matrix elements of the 
Keldysh Green function at different times and its equal- 
time counterparts. To derive the matrix form of the 
GKBA given in p.24p formally and to identify the terms 
which are neglected if one uses this ansatz, we follow 
the derivation by Lipavsky et al.A^ Without any loss of 
generality we will neglect the momentum labels for a mo- 
ment and concentrate on the time dependence only. In 
addition we use the short- hand notation Gu' = [G]tt' for 
the matrix elements in the time labels. We introduce 



KR 



KA 



G^^it,t') = &it-t')G'^{t,t'), (Bla) 
G^^(t, t') = eit' - t)G^{t, t'), (Bib) 



[G"%,,^^GiiZGt,,+eit'^t) 

ft' nt 

X 



dh I dh [Gf,X.t^ + G^tA] G^,. 

J — OO 

(B8) 

Adding Eqs. (|B7P and (|B8[) we obtain the following exact 
integral equation for the Keldysh component of the Green 
function, 

[G%t'^-^ [G^,,ZG^,,-Gf,ZGtt,] 

+ Q{t t') fdh f dh G« [S^,,G^,, + ^iUGtt'] 

Jt' J~oo 

+ Qit' ~ t) fdh f dh [GfuKt, + G^J,S,^tJ Gt'- 

Jt J -oo 

(B9) 
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To rewrite this equation in a more compact form we in- 
troduce the functions 

/oo 
dt,{^l_GfJ + J:f,^Gt,), (BlOa) 
-OO 

/oo 
dh{Gf,^^j:t, + GlJ:^^,,). (Biob) 
-oo 

Then we may write 



[G^%,, = ~^G^,,ZGi^,, + e{t t') I dt,G?,^W^^,, 



-iG?„ZGf^„-' ' M^nR^^rR 



dtiGf,,Wll^, 



dtiGu^ZF^^, 
'[G''ZF%t>, 



[G 



KA-i 



(Ell) 

iGii ZGi, + Q{t' - t) I dt,W,i Gl,. 



iGf,ZGi, 



dhW.iG^,, 



dtiFfiZGf., 



W'^ZG^] 



and hence 



G 



K 



-G^ZF^ + F^ZG^ 



(B12) 
(B13) 



Here the retarded and advanced component of the distri- 
bution function matrix is defined by 



F^ = F° + ZW^, 



(B14a) 
(B14b) 



with the time-diagonal part given by 

[F%t,=i5{t-t')G''{t,t). 



(B15) 



One easily verifies that the blocks have the following sym- 
metries, 



(B16a) 
(B16b) 
(B16c) 
(B16d) 











(^pRf 









The above relations are all exact. Comparing Eq. (j2.15p 
with Eq. (|B13p . we conclude that the parametrization in 
Eq. (|2.15p is indeed correct, and that 



F = ZF"Z. 



(B17) 



The GKBA amounts to retaining only the diagonal part 
F^ of the distribution function. Then the matrix ele- 
ments of the general relation (|B13p reduce to 



G^{t, t') = -i[G"it, t')ZG^{t\ t') - G^{t, t)ZG^{t, t')]. 

(B18) 

This is identical to the GKBA ansatz (|5.2ip which was 
used to study the toy model. Repeating the above calcu- 
lation in the same fashion including the full momentum 
dependence leads to the relation p.24p . 



APPENDIX C: SECOND ORDER SELF-ENERGY OF THE TOY MODEL 



In this appendix we explicitly give the matrix elements of the non-equilibrium self-energies S2(i, t') of our toy model 
introduced in Sec. |V]to second order in the interaction. Ignoring Hartree type of diagrams which arc implicitly taken 
into account by imposing self-consistency in the first order calculation, the non-equilibrium self-energy to second order 
in the interaction is in the contour basis G {+, — }) given by 



o 2 / 

-Zu pp 



npp' fr<pp'\2 1 'yr'PP' r<pp' cPP' 


f^Pp' ((^PP'\2 


I <)ripp' ripp' i-^pp' 
r aa ^ aa ^ aa 


ripp' (r'PP'\'^ -L onpp' cpp' cpp' 




I nr<pp' npp' npp' 



(CI) 
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where the time labels of all Green functions are {t,t'). Using the relations (|2.34af[2.34cl) we obtain for the normal 
part of the matrix elements of the self-energy in the Keldysh (RAK)-basis, 

2 

(C2a) 

E^j,a(i, t') = {replace A'ln the above expression for E|^^g(i, t') } , (C2b) 

2 

5^2^aa(^)*') = ^{G^[(Gajj)^ -I- (Gg^)^ -I- (G^)^] + 2[G^gGa^ -I- G;faGa^]G|^ 

1 nr^K \f^R _i_ _L 1 j_ OC^ G^ + (^^^ 1 4- 9G^ \G^ G^ 4- G^ 1 I 

"I" ^aal^aa^aa ' ^aa^dd ' ^ aa^ ddi ' ^ddl^dd^dd ' ^aa^ddi ' ^ddl^dd^dd ' ^dd^ddi [ 
2 

2 "1^ aa^v aa-' v^da) J ' ^aa^aa^aa ^ '^^daV^aa^dd ^ aa^ dd\ ~^ ^daV^aa^dd ^ ^ aa^ dd\ j 7 

(C2c) 

where Gl.^,{t,t') is the matrix element of the matrix G^ defined in Eq. p.l6p . i.e., 

Gi,, (t, t') = i[Gt' (t, t') - Gt. (t, t')] , (C3) 

and we have used the fact that G^{t,t')G^{t,t') = 0. The corresponding self energies ^^^^(i,^'), E^g^(i,t'), and 
^faaC^' i') can be obtained by simply exchanging a O a in the above expressions. The anomalous components of the 
self-energy are 

2 

'^2,aa{t^t') = ^|Gfa[(Gfg)^ -|- (G^)^] -l- 2G^GaaG^ -|- '2G^g[G^g^G^^ + G^d^da] + 2G^[GfaG^ -I- G^Gg^jj, 

(C4a) 

^2^aa(^i = {replace R^ Am the above expression for E|^^^(i, t')} , (C4b) 

2 

^2^aa(^:*') = ^|G^[(Gaa)^ + {G^df + (Ggg)"] + 2[G^^G?;^ + G^^GaaJGl^ 

+ ^'-^aal'-^aa'-'aa + '-faa'-^aa + '-'aa'-^aaJ + ^"-^aa ['-^aa'-^aa + '-'aa'-^aaJ + ^'-^aal'-^aa'-'aa + '-^aa'-^aaJ J 

2 

2 aal-v aa/ v^dd) J ' ^aa^aa^aa ^ ~^^aaL^aa^aa ^aa^aaJ ^ ^aaL^aa^aa ^aa^aaJ j • 

(C4c) 

Finally, the conjugate anomalous self-energies Ef-aa(t,t'), E^gg(t,t'), and E2^gg(t,i') can be obtained by exchanging 
a O a on both sides of Eqs. (|C4allC4c| . In order to calculate the "out-scattering term" ()3.20bp in the kinetic equation, 
we need only the difference between retarded and advanced self-energies, which to second order in the interaction can 
be written as 

2 

^2,aa(^i*') = *p2!aa(^i*') " ^2^aa(*j^')] = | G^g [(Gg^j)^ — (Gg^)^] - 2G^Ga^Ga^ 

+2Ga(i[G^^Ggg — G^G^] — 2G^[G^^G^ -I- G^Ggg]|, (C5) 

2 

^2,aa(*'^') = *p2^aa(*'*') ~ ^^aa(^:^')] — | G^^ [(Ggg)^ — (Ggg)^] - 2G^GggGga 

+2Gaa[G^aGaQ — G^gG^^] — 2G|^[GQaG^ + G^Gg^jj. (C6) 

The functions E2 aa(^' ^-^id E2 aa(t, t') can again be obtained by exchanging a -H- a in the above expressions. 
APPENDIX D:' SOLUTION OF THE TWO-TIME Markov approximation, and the neglected renormaliza- 
KINETIC EQUATIONS FOR THE TOY MODEL tion of the retarded and advanced propagators. For sim- 
plicity, we focus on the kinetic equations for our toy 
model with the perturbative second-order self-energy. 
In this appendix we exami ne th e validity of the three We compare different combinations of these approxima- 
approximations made in Sec.|VB|to obtain a closed sys- tions and their infiuence on the rcsuhs. For our toy model 
tem of equations for the equal-time Keldysh-Green func- we can obtain the quantum dynamics without relying 
tion of our toy model: the Kadanoff-Baym ansatz, the 
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FIG. 8. (Color online) Time evolution of the diagonal dis- 
tribution function for the second-order self-energy without 
Markov approximation. We compare the result with GKBA 
(top) and the full two-time result without GKBA (bottom). 
For the retarded and advanced propagators we used the renor- 
malized (solid line) or the free propagators (dashed line). The 
parameters and initial conditions are the same as in the mid- 
dle panel in Fig.[T] The dashed line is again the exact solution. 



on any of these approximations by solving the two-time 
quantum dynamic partial differential equations numer- 
ically. Technically, this is almost as simple as solving 
a system of ordinary differential equations in the equal- 
time formalism, except that now we have to propagate 
in two different time directions t and t' . The collision in- 
tegrals were calculated numerically using the trapezoidal 
rule. To shorten the presentation, we will concentrate on 
the dynamics of the normal pair-correlator with the in- 
teraction strength u/e = 0.1 corresponding to the middle 
panel in Fig.|31 The pumping strength I7I and the initial 
conditions are the same as before. 



To begin with, we have solved the kinetic equation 
(|O0l) for using the GKBA but without Markov ap- 
proximation. The results are shown in the upper panel 
of Fig. |S] where we compare two different variants, de- 
pending how the retarded and advanced Green functions 
entering the collision integral in Eq. (|5.20p are handled. 
In the first case we solved the kinetic equation ()5.20|) 
for together with the equations for the renormalized 
retarded and the advanced Green functions, which can 
be obtained from Eq. p.lO|) by simply omitting the mo- 
mentum labels. The set of kinetic equations was solved 
self-consistently using only the GKBA in the collision in- 
tegrals. This combination turned out to be quite unstable 
and the solutions of the equations diverge slightly above 
the time /it/7r = 2. We have checked that the divergence 
is not an artifact of our grid discretization. In the sec- 
ond case wc used the free retarded and advanced Green 
functions in the collision integral of the kinetic equation 
(|5.20p . The solution is again stable but shows an irregular 
dynamics in comparison with the exact solution. Com- 
pared to the analogous results relying in addition on the 
Markov approximation shown in Fig. [3] we did not find 
any improvement. 

Next, we additionally avoided the GKBA. In the lower 
panel of Fig. [8] we compare the two different variants 
using the full renormalized or the free retarded and ad- 
vanced propagators in the collision integral. In the first 
case with the full renormalized quantities, the set of ki- 
netic equations p.lOI3.11|3.12p without momentum la- 
bels was solved simultaneously. The dynamics is sta- 
ble, but the oscillation amplitude disappears nearly com- 
pletely. In the variant with the free advanced and re- 
tarded Green functions the dynamics changes completely 
and the solution shows large oscillations of the pair- 
correlator amplitude. Again we do not observe any im- 
provements towards the correct solution. 

In summary, despite the known limitations of the 
GKBA and the Markov approximations, we have not 
found any improvements in the full two-time approach. 
Note that a similar comparison of the GKBA with full 
and free propagators, respectively, was performed for 
semiconductors in Ref. [1^ , with the conclusion that the 
full GKBA performs rather well. At this point it is not 
clear to us if the different behavior in the present paper 
is due to a breakdown of standard perturbation theory or 
simply a special feature of the toy model, which contains 
no intrinsic dissipation. 
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